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1.  INTRODUCTION 


Currently,  a  number  of  diverse  efforts  are  underway  toward  modeling  and  diagnostics  of 
the  electrothermal-chemical  (ETC)  gun.  These  efforts  have  been  initiated  primarily  in  the  past 
two  years,  include  Government  (Army,  Navy,  DNA,  and  DOE),  university,  and  industry,  and 
are  funded  by  both  private  and  Government  sectors. 

The  three  (Army,  Navy,  and  DNA)  major  Government  programs  associated  with 
development  of  ETC  technology  have  target  dates  of  FY92  for  assessment.  Thus,  a  need 
exists  to  increase  and  encourage  progress  toward  understanding  the  dominant  physical 
mechanisms  in  the  ETC  gun,  which  hopefully,  will  result  in  improved  control  of  the  interior 
ballistic  process. 

As  a  means  of  addressing  the  above  concerns,  a  JANNAF  workshop  on  Electrothermal- 
Chemical  Modeling  and  Diagnostics  was  held  July  9-11,  1991,  at  the  U.S.  Army  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  MD.  the  objectives  of  the  workshop  were  to 
assemble  experts,  drawn  from  gun,  plasma  physics,  engineering,  and  related  disciplines  from 
Government,  industry,  and  academia  to  examine  theoretical  methodologies  and  experimental 
approaches  and  data,  and  to  review  and  evaluate  the  present  state-of-knowledge  In  the  ETC 
gun.  Specifically,  the  workshop  objectives  were  to: 

-  Survey  methods  of  modeling  interior  ballistic  process,  particularly  the  interaction  of  the 
plasma  and  the  work  fluid. 

-  Summarize  the  areas  of  agreement  and  determine  diagnostic  experiments  needed  to 
validate  hypotheses  and  provide  input  for  models. 

-  Identify  diagnostic  experiments  which  may  impact  modeling. 

-  Assess  the  state  of  plasma  modeling  and  diagnostics  for  ETC  guns. 

-  Identify  gaps  in  experimental  and  theoretical  investigations. 

-  Recommend  future  ETC  gun  research  areas. 
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Workshop  participants  jointly  summarized  current  modeling  and  diagnostic  efforts  in  the 
U.S.  and  experimental  measurements  needed  to  improve  ETC  models.  Their  summary  is  in 
the  form  of  the  following  charts  (see  vol.  1): 

a.  Current  ETC  Modeling  Activities  in  the  U.S. 

b.  Diagnostic  Measurements  Desired  by  ETC  Modelers 

c.  Current  ETC  Diagnostic  Activities  in  the  U.S. 

d.  Use  of  Diagnostic  Measurements  Desired  by  ETC  Modelers 

The  dialogue  between  modelers  and  experimentalists  will,  hopefully,  provide  a  common  focus 
for  future  work. 
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ELECTROTHERHAL- CHEMICAL  GUN  PROGRAM 


D.  A.  Benson,  P.  Cahill,  S.  N.  Kempka,  and  R.  L.  Woodfin 
Sandla  National  Laboratories 
Albuquerque,  NM  87185-5800 


ABSTRACT 


The  Sandia  National  Laboratories  ETC  Gun  Program  seeks  to  achieve  a  better 
understanding  of  the  fundamental  processes  at  work  in  such  a  device.  An 
experimental  apparatus  resembling  a  gun  chamber  with  an  external  "capillary"  for 
the  generation  of  a  plasma  which  mixes  with  a  working  fluid  in  the  chamber  has 
been  constructed.  A  series  of  detailed  models  of  the  process  is  being  developed. 
Several  diagnostics  techniques  are  being  developed  to  better  measure  the  parameters 
of  the  process.  Several  possible  fuel  materials,  are  being  investigated.  Future 
plans  include  a  more  "gun- like"  apparatus  with  variable  geometries. 
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Critical  needs: 

*New  working  fluids 
*Well  instrumented  firings 
*Plasma  /  working  fluid  diagnostics 
Bore  temperature  and  erosion  studies 
Projections  of  growth  of  power  supplies 
Major  unknowns: 

*Physics  &  chemistry  of  mixing 
*Gas  generation  rate  t/s  electrical  input 
*Gun  thermal  environment 
Heat  transfer  &  erosion 
Muzzle  reactions  &  weapon  effluents 
Ballistic  temperature  coefficient 
‘Reproducibility 


Sandia  National  Laboratories 
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Physics  Of  plasma  arc  sources 

Physics  &  chemistry  of  mixing  processes 

Modeling  of  the  electric  /  thermal  /  chemical  process 

Properties  of  propellant  materials 

Gun  chamber  diagnostics 

Improvements  in  plasma  density  /  mixing 
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Integrated  experimental  &  analytical  effort 


Produces  &  verifies  detailed  physical  models 


Supports  Army  modeling  of  interior  ballistics 


Develops  new  diagnostics  for  mixing  &  combustion 
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Goal  is  to  be  able  to  identify: 

plasma  penetration  deptfi  &  volume 
& 

location  of  species 

under  gun  conditions 
as  a  function  of  time. 


Diagnostics  will  measure: 

Plasma  development 

Temperature 

Species 
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Experimental  chamber 
To  100,000  psi 
40KJ 

Heavily  instrumented 
Contained  products 
Ablative  capillary  liner 
Analytical  model 

2-D,  axisymmetric,  compressible  flow 
Dissociated  &  ionized  species 
Electrical  conductivity 
Gas  equation  of  state 

Radiative  energy  transport  by  diffusion  approximation 


Sandia  National  Laboratories 


Plasma  /  working  fluid  interaction 

Flow 

Mixing 

Areas  of  reactivity 
Entrainment  processes 

Ablation  of  capillary  material 

Mass  /  momentum  addition  to  plasma 
Reactivity  alteration  in  accumulator 

Constitution  of  plasma 

Mixture  of  ablative  material  &  working  fluid 


RLW  7/3/91 


Sandia  National  Laboratories 


mwm 


Electrothermal  Injector  Studies 


Develop  ETC  technology  base'for  future  systerrt.assessment 
H  Fundamental  electrothermal  chemical  measurements  . 

H  Materials  developrnent  ,  .  .  . 

■  Modeling 

H  Coordinate  with  DoD  and  industry  efforts 


simplified  bUN  DESIGN 

REAR  capillary  ■■  OXIDIZER 

ELECTRODE  1^  LINER  (FUEL)  FLUID  PROJECTILE 


Objectives 

Hf:  Develop  closed  vessel  electrothermal  test 
0  Acquire  data  for  capillary  model  development 
H;  Improve  efficiency  and  controllability  of  ETC  devices 
H  Enhance  capillary  .mass  transfer 


ELECTROTHERMAL  INJECTOR  TEST 
CAPILLARY  TEST 


^Evaluate  reactive  materials  for  use  in  capillary 
H  Study  capillary  turbulence,  heat  transfer,  erosion 
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Electrothermal  Injector  Experimental j^Xpparatus 


Capacitor  Discharge  System 
40  kJ,  250  kA,  10  kV 


■  cleclrlcal 

pfoba 


pressure  (Ittlnge 

accumulator 
/  chamber 


H  Liner  Materials 
Polyethylene 
Delrin 

Nitrated  Polyethylene 

SI  Diagnostics  " 

Arc  Impedance 
Electrical  Power  Input 
Accumulator  Chamber  Pressure 
Capillary  Pressure  vs  Time 
Total  Mass  Rernoval  from  Capillary  Liner 

H  Post-Test  Particulate  Analysis 
Scanning  Electron  Microscopy 
Energy  Dispersive  X-Ray  Analysis 


capillary 


'  pressure 
transducers 


Electrothermal  Injector  Cross-Section  View 

H  Experirriental  Setup 

Accumulator  Volume  14  cm^ 
r 

Arc  Length  4.3  cm 
Capillary  Inside  Dia.  0.58  cm 
Max.  Pressure  100  kpsi 
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mixing 
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Accomplishments  on  Electrothermal  Injector  Project 


■  Developed  closed  vessel  test  system  for 

the  Dynamic  Heating  Facility 

■  Conducted  initial  tests  with  polyethylene 

(10  to  80  mg  ablated) 

H  Calculated  optical  absorption  to  predict 
radiative  paths  under  ET  conditions 

H  Performed  hydrodynamic  calculations  to 
show  characteristic  times  for  capillary 

H  Prepared  and  characterized  nitrated 

polyethylene  materials  (PEVN)  for  liners 

H  Demonstrated  enhanced  arc  induced 
ablation  with  PEVN  liners 


arc  power  dissipated 


0.0  0.2  0.4  0.6  0.8  1.0 

Time  (ms) 


accumulator  chamber  pressure 


Laboratories 


Model  of  injector  /  accumulator  | 


Electric  discharge 
Plasma  generation 
2-D  fluid  flow 

Adaptation  of  existing  codes 
Mixing,  combustion,  heat  transfer 
Geometric  sensitivity  studies 

Model  of  ablation  process 

Code  models  measurements 
Input  to  2-D  model 
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1-D  Radial  Model  of  an  ETC  Capillary 


Mivina  Ph3m>v».r 


•Objectives 

1 .  Examine  extent  of  radial  variations  in  plasma 

-  radial  variations  influence  mixing,  instabilities 

2.  Validate  plasma  transport  properties  for  use  in  2-D  hydro  model 

-  diffusive  radiation  transport  (Rosseland  mean  free  path) 

-  electrical  conductivity  (electron-ipn,  electron-neutral  collisions) 
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Energy  Balance 


hj 


Ohmic  Heating 
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Assumptions  in  1-D  Radial  Model 

•  Steady-State,  LTE,  (collision  rate)'''  « 


’  Axially  uniform  velocity  and  temperature 

•  Axial  flow:  Exit  capillary  at  local  sonic  velocity 

•  Radial  flow:  obtained  from  conservation  Of  mass 


•  Radial  pressure  gradient  from  radial  momentum  equation  (mhd  included) 


•  Uniform  Eiectrical  field 


•  Saha  equations:  C,  C'^,  H, 
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Computed  Radial  Distributions 


Find  solution  that  matches  SNL  experiniental  data 
-  Peak  values:  4d  'kAmps,  100  MW,  000  V/cm 
temperature  varies  from  4.8  eV  to  2.5  e-V  across  capillary  radius  (0.25  cm) 

5  ^ — I — ' — ' — ■ — ' — I — ' — I — --“I — ]■ '  ■■  — I — ^ ^ 

[  i  •  Species  also  vary 

‘f  ■;  significantly  with  radius 


V  V  V \'VV] 


Centerline 


0,10  0.15 

Rodius'  (cm) 


Uniform  temperature 
-  Power  =  oT^*  2tiRL 


Capillary;,Raclius 


RLW  7/3/91 


Sandla  Matronal  Laboratories 


1-D  Radial  Moddl  Conclusions 


•  Radial  model  provides  insight  into  capllary  physics: 

-  Radial  variations  in  plasma  tempelature,  species  are  large 

-  Radial  variations  effect  mixing  instabilities,  control  issue 

•  Modeling  issues  to  be  resolved: 

-  Ionization  energy  depends  on  density  (constant  value  used) 

-  Electrical  conductivity  of  plasma  for  high  densities 


Sandia  National  Laboratories 


2-D  Hydrodynamic  Model 


•  Include  capillary  model  in  2-D  hydro  model  to  examine  scaling,  mixing 

•  Time-step  issue 

-  Explicit:  expensive,  accurate 

-  Implicit:  cheaper,  large  numerical  diffusion,  dispersion 

•  Explicit:  time-split  diffusive  and  convectivp  transport:  different  time  scales 

•  Flux  Corrected  Transport  for  convective  transport 

•  Test:  Numerical  Simulation  of  Convection:  Propagating  Wave 
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•  Ret:  Oran,  E.  S.  and  J.  P.  Boris,  •‘Detailed  Modeling  of  Combostion  Systems,"  Prog. 
Energy  Combust  Sd..  1981 
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Reactive  injector  liners 
Graded  liners  /  current  profiling 
Ablation  rate  measurements 
Modeling  of  ablation 


Observations  of  density  vs  mixing 


Developing  copolymer  reactive  liners 

.  Nitrated  ethylene  /  PVA  copolyn|er 
For  increasing  plasma  density 
For  use  with  HAN 

Ablation  rate  calibration 

From  shock  tube  experiments 

Gun  simulator  testing  of  above  materials 

For  plasma  density  effects 
For  mixing  improvements 
For  energetics  /  kinetics 
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Ablation  Data  from  Electrothermal  Infector  Tests  for  Three  Materials 


Test  Energy  (kJ) 


B  Washer  tests  (PEVN),  preliminary  safety  testing 

B  Full  liner  tests,  normal  testing  configuration 

B  Specific  Energy  —  Elec.  Energy  to  arc  /  Total  liner^ass  removed 
(Small  value  — >  Large  mass  removed) 
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Studying  fundamental  processes 
Provide  inputs  to  higher  level  models 
Supports  BRL  efforts 
Will  answer  some  scaling  issues 
Methodical,  low  risk  approach 


DIAGNOSTICS  DEVELOPMENT  FOR  THE  ETC  PROGRAM 


Donald  W.  Sweeney,  Steven  Vosen,  Jeffrey  Gray,  and  Robert  Armstrong 

Sandia  National  Laboratories 
Livermore,  CA  94551-0969 


ABSTRACT 

Control  of  the  ET  process  is  governed  by  the  proper  coupling  of  plasma 
energy  to  chemical  energy.  Optimization  requires  that  the  mixing  and  combustion 
process  be  understood.  The  physics  and  chemistry  of  the  mixing  of  a  plasma 
injected  into  a  confined  liquid  will  be  studied  by  combining  the  techniques  of 
high-speed  cinematography  and  emission  spectroscopy.  Variation  of  the  mixing 
chamber  geometry  and  of  the  liquid,  with  increasing  energetic  content  (e.g., 
from  water  to  methanol  to  hydroxyl  ammonium  nitrate) ,  will  allow  the  effects  of 
chemistry  on  plume  development  and  mixing  to  be  studied. 
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DIAGNOSTIC  DEVELOPMENT  FOR  THE  ETC 

PROGRAM 


JANNAF  Workshop 
ETC  Modeling  &  Diagnostics 

O' 


July  9-11,  1991 


Donald  Sweeney,  Steve  Vosen, 
Jeff  Gray  and  Rob  Armstrong 


Combustion  Research  Facility 
Sandia  National  Laboratories 


Will  Review  Three  Components  of  the  CRF  Program 

1 .  Diagnostics  of  plasma  generation  and  fluid  mixing 

2.  Chemical  and  thermal  ablation  studies 

3.  Laser-induced  plasmas  for  ETC  development 


Diagnostics  for  measuring  gas  generation  and  plasma 
behavior  will  be  developed  for  gun  conditions _ 


Development  of  diagnostics  for  use  at  gun 
conditions  is  underway  to  quantify: 

1)  plasma  behavior,  and 

2)  gas  generation  rates. 

Development  will  progress  through  three  pressure 
regimes: 

1)  Low  Pressure  :  P  <  1  kpsi  (FY  91), 

2)  Medium  Pressure:  P  =  1  - 15  kpsi  (FY  92),  and 

3)  High  Pressure:  P  =  10  - 100  kpsi  (FY  93). 


At  low  pressures  we  will  observe  fluid  mixing 
over  a  range  of  pressures 


Pressures:  15-1000 


Process/Location 

Diagnostic  Technique 
•  Measurement 

Gas  Generation/ 
Mixing  Chamber 

Photography  and  image  processing 
•  Deveiop  image  processing  techniques, 
simiiar  to  those  deveioped  for  iiquid 
propeiiants,  to  quantify  images  of  fuii 
fiow-fieid  mixing 

Vary  energy  and  mixing  fiuid. 

*  Determine  effect  of  energy  and  fiuid  on 
evoiution  of  mixing  region. 

Piasma  Formation/ 
Piasma  Cartridge 

n.a. 

(expiosive  charge) 

Pressure  vessel 


Pressure:  0  -  5000  psig 
Volume:  13  liters 
Windows:  4, 100  mm  dia. 


Mixingfimage  processing  experiments  will 
be  done  in  a  low  pressure  apparatus 


Pressure  (psi) 


Pressure  cartridges  release  200  Joules  ini  ~  2  milliseconds 


Time  (ms) 


Relevant  gas  generation  mechanisms  and  plasma 
behavior  will  he  studied  at  medium  pressure 

Pressures:  1-15  kpsi  | 


Process/Location 

Diagnostic  Technique 
•  Measurement 

Gas  Generation/ 
Mixing  Chamber 

Photography/image  processing 
•  Characterize  mixing  over  a 
range  of  conditions. 

(Maximum  pressure  overlaps  with  high 
pressure  studies) 

Plasma  Formation/ 
Plasma  Cartridge 

Time  resolved  emission 

•  Plasma  stability 

Spectrally  resolved  emission 

•  Plasma  temperature 

•  Electron  density 

•  Sheath  temperature 

•  Identify  Species 

Heat  transfer  probes 

•  Radiative  heat  transfer 

•  Conductive  heat  transfer 

A  Medium  Pressure  Apparatus  will  be  used  to  study  plasma 
characteristics  and  flow-field  mixing 


Schematic  of  medium 
pressure  apparatus 


The  transparent  chamber  with  corrective  optics  will  allow 
for  Schlieren  photographic  measurements 


■1^ 


bH 

f/l^ 

WPM^f^^iAijML 

UJU^ 

Signal  (a.u.) 


Fiber  optics  will  provide  spectral  and! or  temporal  information 


The  shock  tube  provides  a  controlled  environment  to 
observe  both  thermal  and  chemical  ablation  processes 


Thermal  boundary 


Full  gas-phase  and 
surface  chemistry 
(CHEMKIN) 


Experimental  Schematic 


Nd ;  YAG-pumped 
Dye  Laser 
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Chemical  ablation  is  initiated  by  seeding 

the  shock  tube 


N2  O  -►  N2  +  o 

The  0-atoms  attack  the  surface  and  generate  OH 

O  +  (-CH2-)„  -►oh  +  (-CHGH2-) 

OH  is  observed  using  LIF 

1“ 


Data  and  model  agree  when  the  model  is  blurred 
to  account  for  instrument  function 


1 


Data  and  model  agree  when  the  model  is  blurred 
to  account  for  instrument  function 


1 


ETC  Combustion  is  being  studied  using 
Laser  -  Plasma  Ablation 


©  Spectroscopy 
©  Imaging 


Camera 


Laser  Pulse 


Laser  -  induced  plasmas  help  discern  various 

ETC  processes 


r  Dense  Plasma 


Thermal  Ablation 


Ablation 


Brittle 

Fragmentation 


Reflected  Shock 


Plasma 


Vaporized 

Material 


1: 

hi 

^  T 

K’lxx.<:i 

Fragments 


Variations  in  Test  Conditions 


8  11  3 

•  Plasma  power  density  ranges  from  10  to  10  W/cm 

•  Samples  of  PE,  PEVN,  etc. 

•  Reflected  shock  conditions  to  >  6000  K,  >  10^  atm 

Kn 

so 

•  Plasma  formation  in  air,  He,  Ar,  vacuum,  and  water 

•  Temporal  resolution  of  20  ns,  Spatial  resolution  of  20 

•  Spectroscopy  from  0.2  to  1.0  ^m 


Laser  -  Ablation  Spectroscopy  Reveals 
Plasma  Species  for  Dijferent  Materials 


A 


740  750  760  770  780 

X  (nm) 


Time  -  Resolved  Imaging  Shows  Ablation  Rates 

for  Different  Materials 


Polyethylene  at  low  pressure 


0  delay 


50  ns  delay 


Ablation  Rates  may  be  related  to  the 
Energetics  of  Materials 


PEVN  at  low  pressure 


0  delay 


50  ns  delay 


Laser  Plasmas  Simulate  Conditions  on  an 
Optically  Accessible  Laboratory  Scale 


o  Highest  power  density  achieved  in 
Sandia  tests  (>  10^®W/cm^) 

o  Pressures  and  temperatures  in  the  same 
range  as  gun  conditions 

©  High  repetition  rate  (10  Hz)  allows  us  to 
perform  a  large  matrix  of  tests  studying 
ablation,  and  possibly  mixing 

o  Provides  a  diagnostics  testbed  for 
full  scale  tests 


Intentionally  left  blank. 
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DEVELOPMENT  OF  AN  UPtfINO/IMPLICIT  COMPUTATIONAL  MODEL 
FOR  THE  ADVANCEMENT  OF  ARMY  ETC  GUNS 

N.  Slnha,  A.  Hosangadl,  and  S.  H.  Dash 
Science  Applications  International  Corporation 
Fort  Washington,  PA 

ABSTRACT 

A  first-principles,  state-of-the-art  computational  model  based  on  the  CRAFT 
Navier-Stokes  solver  is  under  development  for  the  U.S.  Amy  to  simulate  fluid 
dynamic  and  electro-themo- chemical  processes  in  ETC  guns.  The  model  will 
incorporate  physics  submodules  developed  from  diagnostics  by  the  team  of  SAIC, 
BRL  and  SNL  scf^entists  as  well  as  national  and  international  efforts  coordinated 
with  the  Army. 

This  presentation  will  provide  an  overview  of  the  1D/2D/3D  finite-volume 
upwind  (Roe/TVD)/lmpllclt  numerics  in  CRAFT  and  its  strongly-coupled  treatment 
of  combustion  chemistry  and  multi-phase  non-equilibrium  processes  (presently 
limited  to  dilute  particulates).  Current  work  focussed  on  the  ETC  gun  problem 
includes  dynamic  grid  capabilities  with  solution  adaptive  features,  and  research 
towards  incorporating  complex  equations  of  state  into  higher-order  Rlemann  based 
solvers  such  as  the  Roe/TVD  scheme  implemented.  A  building-block  approach  upgrade 
of  CRAFT  to  simulate  ETC  processes  will  be  described  along  with  representative 
calculations  which  exhibit  present  capabilities. 
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DEVELOPMENT  OF  AN  UPWIND/IMPLICIT 
COMPUTATIONAL  MODEL  FOR  THE 
ADVANCEMENT  OF  ARMY  ETC  GUNS 

N.  Sinha,  A.  Hosangadi,  and  S.  M.  Dash 
Science  Applications  International  Corporation 
Fluid  and  Propulsion  Sciences  Operation 
501  Office  Center  Drive,  Suite  420 
Fort  Washington,  PA  19034-3211 
Phone:  215/542-1200  FAX:  215/542-8076 


JULY  9-11,  1991 


Presented  at: 

JANNAF  WORKSHOP  ON  ETC  MODELING  AND  DIAGNOSTICS 

U.S.  ARMY  BALLISTICS  RESEARCH  LABORATORY 
ABERDEEN  PROVING  GROUND.  MARYLAND 


'•  An  Employee-Owned  Company 

SclenoeAflpScaihnshternationalCotpofation 


TOPICS 


•  APPROACH  TO  BE  TAKEN 

•  CRAFT  1D/2D/3D  UPWIND/IMPLICIT  FINITE-VOLUME  CODE 

-  HISTORY 

-  EQUATIONS/NUMERiCS 

-  STRONGLY-COUPLED  CHEMISTRY 

-  COUPLED  ANALYSIS  OF  DILUTE  PARTICULATES 

-  DYNAMIC  GRID 

O' 

•  RELATED  APPLICATIONS 

•  CRAFT/ETC  DEVELOPMENT 

-  ADAPTION  TO  ETC  ENVIRONMENT 

-  COMPLEX  EQUATIONS  OF  STATE/RIEMANN  PROBLEM 

-  GAS/LIQUID  PROBLEMS 


•  PRELIMINARY  APPLICATIONS  (BRL  PROGRAM  INITIATED  6/1/91) 


An  Employoe-Owned  Company 

SckncBAfipScationsIntemiantv^Conxjmtion 


BRL  APPROACH  FOR  ETC  GUN  SIMULATION 


EMPHASIS  ON  UNDERSTANDING  "PHENOMENOLOGY" 
AND  INFLUENCE  OF  THERMOCHEMISTRY  AND  PHYSICS 
ON  GUN  PERFORMANCE 


EXPERIMENTAL  DATA  USED  TO  SUPPORT  DEVELOPMENT 
OF  "FIRST-PRINCIPLES"  CFD  SIMULATION  CODE  AND  TO 
ENHANCE  UNDERSTANDING  OF  PHENOMENOLOGY 

-  NO  ATTEMPTS  TO  "TUNE"  THERMOCHEMICAL/ 
PHYSICAL  PARAMETERS  TO  MATCH  DATA 


BUILDING  BLOCK  APPROACH  FOR  UNDERSTANDING 
PHENOMENOLOGY 


EXPLORE  INDIVIDUAL  EFFECTS 


Army  ETC  Modeling  &  Diagnostics 


mod/d1ag.may91 


FEATURES  OF  CRAFT  FOR  ETC  GUN  SIMULATION 


FEATURE 

1  ADVANTAGE 

1D/2D/AXI/3D  RUN 
OPTIONS 

CAN  DO  FAST  PRELIMINARY  DESIGN  STUDIES  (ID)  AND  DETAILED 
SIMULATION  (2D/AXI/3D)  IN  SINGLE  CODE 

FINITE-VOLUME 

IMPLICIT 

CONSERVATIVE 

FORMULATION 

CONSERVATIVE/FINITE-VOLUME  NUMERICS  YIELDS  ACCURATE  FLUX 
BALANCES  -  EQUATIONS  STRONGLY  COUPLED;  IMPLICIT  NUMERICS 
ELIMINATES  STABILITY/  STEP-SIZE  PROBLEMS  AND  PERMITS 

RESOLVING  NEAR-WALL  BOUNDARY  LAYERS 

ROE/TVD  UPWIND 
NUMERICS 

ROE  IS  REIMANN  BASED  DIFFERENCING  FOR  ACCURATE 
CONVECTIVE/WAVE  TRACKING;  TVD  IS  HIGHER-ORDER  EXTENSION;  NO 
ARTIFICIAL  DISSIPATION  REQUIRED  FOR  STABILITY 

BLOCK  MATRIX 
INVERSION  WITH  INNER 
NEWTON  ITERATION 

ROBUST  INVERSION  PROCEDURE  PERMITS  TAKING  LARGE  TIME-STEPS; 
INNER  ITERATIONS  PERMITS  ELIMINATING  ALL 
FACTORIZATION/LINEARIZATION  ERRORS  WITHIN  STEP 

STRONGLY-COUPLED, 

FULLY-IMPLICIT 

CHEMISTRY 

CHEMICAL  SPECIES  ARE  PART  OF  COUPLED  EQUATION  SET  FOR 
ACCURACY/ROBUSTNESS.  CHEMICAL  SOURCE  TERMS  TREATED  FULLY- 
IMPLICITLY  TO  ELIMINATE  STIFFNESS 

STRONGLY-COUPLED, 
FULLY-IMPLICIT 
TURBULENCE  MODELS 

TURBULENCE  MODEL  EQUATIONS  ARE  PART  OF  COUPLED  EQUATION 

SET  FOR  ACCURACY/ROBUSTNESS;  TURBULENCE  SOURCE  TERMS 
TREATED  FULLY-IMPLICITLY  TO  ELIMINATE  STIFFNESS 

IMPLICIT/UPWIND 
DISPERSED  PHASE 
SOLVER 

PARTICULATE  EQUATIONS  SOLVER  ON  SAME  GRID  AS  FLUID  PHASE 
USING  SAME  IMPLICIT/UPWIND  NUMERICS;  STRONG  FLUID/PARTICLE 
COUPLING  AND  IMPLICIT  SOURCE  TERM  TREATMENT 

DYNAMIC  GRID, 
SOLUTION  ADAPTIVE 
MODULE 

- 1 

GRID  MOVES/DEFORMS  IN  ACCORDANCE  WITH  BOUNDARY  MOTION 
(PROJECTILE)  AND  REDISTRIBUTES  USING  NASA/AMES  SAGE  MODULE  | 

EVOLUTION  OF  METHODOLOGY  IN  CRAFT  CODE 


UPWIND/IMPLICIT  FINITE-VOLUME  METHODOLOGY  (ROE/TVD)  APPLIED  TO  EXTERNAL 
AERODYNAMIC  FLOWFIELDS  ~  1986-1988.  Reference:  Walters  and  Thomas,  "Advances  in  Upwind 
Relaxation  Methods,"  State  of  the  Art  Surveys  on  Computational  Fluid  Mechanics.  ASME  Pub.,  1988. 

EXTENSIONS  OF  ROE/TVD  METHODOLOGY  TO  MULTI-COMPONENT/REAL  GAS  FLOWS  WITH 
CHEMICAL/THERMAL  NONEQUILIBRIUM  ~  1988-1989.  Reference:  Liu  and  Vinokur,  "Upwind 
Algorithms  for  General  Thermo-Chemical  Nonequilibrium  Flows,"  AIAA  Paper  89-0201,  Jan.  1989. 

EXPLORATORY  STUDIES  OF  ROE/TVD  NUMERICS  FOR  UNSTEADY  FLOWS,  ACCURATE  WAVE 
PROPAGATION  IN  DUCT  (RESONANCE  TUBE)  ~  1989.  Reference:  Ridder  and  Beddini,  "Time- 
Accurate  Finite-Volume  Method  for  Propulsion  Chamber  Flows,"  AIAA  Paper  89-2SS4,  July  1989. 

EXTENSIONS  OF  ROE/TVD  METHODOLOGY  FOR  PROPULSIVE  FLOWS  WITH  COMBUSTION 
CHEMISTRY,  ADVANCED  TURBULENCE  MODELS  ~  1989-1991.  Reference:  Dash,  "Advanced 
Computational  Models  for  Analyzing  High-Speed  Propulsive  Flowfields,"  1990  lANNAF  Propulsion  Meeting. 
Oct.  1990. 

EXTENSIONS  OF  ROE/TVD  METHODOLOGY  FOR  UNSTEADY  ROCKET  PLUME  PROBLEMS  WITH 
MULTI-PHASE  NONEQUILIBRIUM  -  1989-1991.  Reference:  Sinha,  Hosangadi  and  Dash,  "The  CRAFT 
NS  Code  and  Preliminary  Applications  to  Unsteady,  Reacting,  Multi-Phase  Plume  Flowfield  Problems,"  1991 
JANNAF  Plume  Meeting.  May  1991. 

EXPLORATORY  STUDIES  OF  ROE/TVD  METHODOLOGY  FOR  GUN  FLOWFIELDS  WITH 
DYNAMIC  GRID  EXTENSIONS  ~  1990-1991.  Reference:  Hosangadi  and  Sinha,  "Development  of  a 
Dynamic  Grid  Navier-Stokes  Solver  for  Computing  Electro-Thermo-Chemical  Gun  Conflgurations," 
SAIC/TR-88,  Dec.  1990. 
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SYNERGISM  OF  CRAFT  WITH  OTHER  DoD  &  NASA  EFFORTS 


NUMERICS  IN  CRAFT  REPRESENTS  CURRENT  STATE-OF-THE-ART  FOR  VISCOUS, 
CHEMICALLY-REACTING  FLOWS 


ALL  CURRENT  NATIONAL  PROGRAMS  ARE  IMPLEMENTING  SUCH  NUMERICS  FOR  REAL 
PROBLEMS 

NATIONAL  AEROSPACE  PLANE  (NASP)  PROGRAM 
(GASP  CODE,  USA  CODE,  SCHAFT/CRAFT  CODES,  ...) 

-  SCRAMJET  PROPULSION 

SDIO  PROGRAM  (GASP  CODE,  SCHAFT/CRAFT  CODES) 

MISSILE  JETS/PLUMES 

HIGH-SPEED  RESEARCH  PROGRAM  (HSRP) 

AIRCRAFT  JETS/PLUMES 


6.1  RESEARCH  PROGRAMS  USING  THIS  CLASS  OF  NUMERICS  AS  TOOL  TO  EXPLORE 
PM /SICS 

COMBUSTION  INSTABILITY  IN  ROCKETS/RAMJETS 
-  DIRECT  OR  LARGE  EDDY  TURBULENCE  SIMULATION 


SAIC  IS  IMPLEMENTING  VERSIONS  OF  CRAFT  TO  SUPPORT  SUCH  PROGRAMS 


s 


E 


1  CODE  NAME 

CRAFT/JR 

CRAFT/GF 

CRAFT/RNP 

CRAFT/LU 

APPLICATION 

HIGH-SPEED  JET 
RESEARCH 

ETC  GUN  FLOWFIELDS 

ROCKET  NOZZLES/ 
PLUMES 

NUMERICAL 

RESEARCH 

SPONSOR 

NASA/LaRC,  HSRP 

BRL 

MICOM 

INTERNAL 

RESEARCH 

EQUATIONS 

1D/2D/AXI/3D 

1D/2D/AXI/3D 

1D/2D/AXI/3D 

1D/2D/AXI/3D 

CHEMISTRY 

PERFECT  GAS,  TWO- 
STREAM 

IMPERFECT  GAS, 
COMBUSTION 
CHEMISTRY,  MULTI¬ 
PHASE  CHEMISTRY 

SPF  FINITE-RATE 
CHEMISTRY,  TWO- 
STREAM,  PERFECT 

GAS 

PERFECT  GAS 

TURBULENCE 

kc.  COMPRESSI¬ 
BILITY  EXTENSIONS, 

ARS  EXTENSIONS 

kc 

kc/Chien, 

COMPRESSIBILITY 

EXTENSIONS 

kc 

PARTICULATES 

NONE 

EQUILIBRATED 

MIXTURE,  NONEQ. 
LIQUID  PROPELLANTS 

NONEQUILIBRIUM 

NONE 

SOLUTION 

IMPLICIT/UPWIND 
(Roe/TVD)  STRONGLY- 
COUPLED  FLUID/ 
SPECIES/  TURBULENCE 
(8x8  BLOCK 
INVERSION), 
TIME-ACCURATE  OR 
TIME-ASYMPTOTIC 

IMPLICIT/UPWIND 
(Roe/TVD)  ON 

DYNAMIC  GRID,  TIME- 
ACCURATE  SOLUTION, 
STRONGLY-COUPLED 
EQUATIONS,  VARIABLE 
MATRIX  SIZE 

IMPLICIT/UPWIND 
(Roe/TVD)  ON  FIXED 
GRID,  TIME- 
ACCURATE  OR  TIME- 
ASYMPTOTIC, 
STRONGLY-COUPLED 
EQUATIONS, 

VARIABLE  MATRIX 

SIZE 

LU  UPGRADE  FOR 
ROBUSTNESS, 
FASTER 
CONVERGENCE 
(CFL  -  25-50) 

NEW  WORK 

EXPLORING  NEW 
TURBULENCE  MODELS. 
NEW  NONREFLECTfVE 

BC 

INCLUSION  OF  LIQUID 
PHASE,  GAS/  LIQUID 
INTERFACE,  DROPLET 
FORMATION  MODEL 

ADAPTIVE  GRIDDING 
FOR  UNSTEADY 

MULTI-  PHASE  FLOWS 

REWRITE  OF  CODE 
STRUCTURE  TO 
OPTIMIZE  LU 
STORAGE  FOR  3D 

Science  Apfjiications  tntematknat  Corpomtion 


DYNAMIC  GRID  COMPUTATIONS 


j  Qdv  -  j  Qdv  +  a  •  f  dsdt  =  j|**J^^^pdvbdt 


Term  I 


Term  U 


Term  in 


Term  I:  PRIMARY  CONSERVED  VARUBLES 

Term  U:  FLUX  TERMS 

Term  lO;  SOURCE  TERMS 

V(t):  TIME  VARYING  VOLUME 


TREAT  GRID  AS  PURELY  GEOMETRIC  QUANTITY,  I.e.,  DEnNITION  OF 
CONTRAVARIANT  VELOCITY  UNCHANGED 

EVALUATE  F1.UX  TERM  BY  ASSUMING  GRID  IS  HELD  CONSTANT  AT  TIME- 
AVERAGED  VALUE 


METHODOLOGY 


EVALUATE  SOURCE  TERM  FOR  THE  GRID  AT  THE  NEW  TIME  LEVEL,  THUS 
GIVING: 


-Q"  Wi  _ 


V"**  +  Q 


+  F“"'SAt  = 


F  =  a*P,  and  S  is  the  time  averaged  metric 


NUMERICS 


SPECIFY  GRID  MOVEMENT  AS  A  FUNCTION  OF  TIME 
EVALUATE  F„  AS  A  THIRD-ORDER  UPWIND  BASE  FLUX 
SOLVE  LINEARIZED  IMPLICIT  OPERATOR  USING  ADI  PROCEDURE 
USE  NEWTON  ITERATION  FOR  HIGHER  ORDER  ACCURACY  IN  TIME 


•  An  Empkv^e-Owned  Compenf 
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SHOCK  TUBE  COMPUTATIONS 


PRESSURE 

10«1  PRESSURE  RnTIO 


OMtaJ)  l£VCtS 
0.079)0 
O.lOQOO 

S:^ 

0.179)0 


S:^ 


J79)a 

o.SGom 

0.939)0 

O.9SQD0 

0.979)0 

o.iomo 

0.139)0 
o.Isooo 
0.479)0 
n  soom 

«l  'i'  i»:l 
Q.CK'UU 

n. eaHM 
a.eam 
0.67910 
0.70X10 
O.73SD0 

o. 7iSo 
0.775)0 
O.6(U)0 

0.679)0 


Science  Applications  International  Corporation 


10^  IS.0  20.D  ^.0  30.0  35.0  40.0  45.0  9.0  S.0 


-  PRESSURE  VARIATIONS 


PRESSURE 

lOOOit  PRESSURE  RnTIO 


p  CONTOUR 


s 


OLUTION  ADAPTIVE  GRID  GENERATION 


•  NASA/AMES  RESEARCH  CENTER  SAGE  2D/3D  CODE 
"A  Simplified  Self-Adaptive  Grid  Method,  SAGE,"  NASA 
TM-102198. 

•  MODIFIED  VERSION  OF  TECHNIQUE  DEVELOPED  BY 
Diewert  and  Nakahashi  of  NASA/AMES  (AIAA  PAPER 
85-1525,1985) 

•  BASED  ON  VARIATIONAL  PRINCIPLES 

•  PROCEDURE  ANALOGOUS  TO  APPLICATION  OF 
TENSION  AND  TORSION  SPRING  FORCES  PROPOR¬ 
TIONAL  TO  LOCAL  FLOW  GRADIENT  AND  FINDING 
EQUILIBRIUM  (REDISTRIBUTED)  POSITION  OF  THE 
RESULTING  SYSTEM  OF  GRID  POINTS 

•  USER-DEFINED  WEIGHTED  COMBINATION  OF 
ADAPTION  VARIABLES 

•  CURRENTLY  SAGE  BEING  OPTIMIZED  FOR  ETC 
PROBLEMS  TO  YIELD  "BEST"  COMBINATION  OF 
ORTHOGONALITY,  SMOOTHNESS  AND  CLUSTERING 


ScmtoBAfipBeatiomintmnationalCofpontion 
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H/JL 


TRANSIENT  UNDEREXPANDED  JET  INJECTION  INTO  MIXING  CHAMBER 

INITIAL  PRESSURE  RATIO  =  5 


0.0  1.0 
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FEATURES  OF  THE  PAkTICLE-CLOUD  SOLVER 


FEATURES 


EQUATIONS 


_ PARTICLE  CLOUD  SOLVER _ _ 

•  1D/2D/3D  TIME  MARCHING  PARTICL&CLOUD  EQUATIONS 

•  FULLY  COUPLED  PARTICLE  DENSITY,  VELOCITY,  AND  TEMPERATURE  EQUATIONS 

•  DRAG  TERMS  MODELED  FOR  NON-EQUILIBRIUM  PARTICLE  INTERACTIONS 

•  DRAG  RELATIONS  OF  HERMSEN  AND  HENDERSON  AVAILABLE 

•  PHASE  CHANGE  OF  PARTICLES  INCORPORATED _ 

•  LIMITING  PARTICLE  STREAMLINE  CAPTURED  AS  A  DISCONTINUITY  IN  THE  PARTICLES 
CLOUD  DENSITY 

•  ROE  TVD  UPWIND  ALGORITHM  FULLY  CONSISTENT  WITH  GAS  PHASE  NUMERICS 


NUMERICS 


•  IMPLICIT  TREATMENT  OF  PARTICLE  SOURCE  AND  FLUX  TERMS 


BOUNDARY 

CONDITIONS 


ADI  BLOCK  INVERSION 

TIME  ACCURATE/TIME  ASYMPTOTIC  NUMERICS 

NEWTON  ITERATION  TO  REMOVE  APPROXIMATE  FACTORIZATION  ERRORS 
WALL  WITH  PARTICLES  STICKING 

WALL  WITH  PARTICLES  SLIPPING 

INFLOW 

OUTFLOW 
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Finite-Duration  Transverse  Particle 
Injection:  Density  Contours  of  Particle 
Packet  at  Four  Time  Intervals 


MULTI-PHASE  PERPENDICULAR  JET  INJECTION 
PARAMETRIC  VARIATION  OF  PARTICLE-SIZE 


GAS  FLOW 


M  =  2 
P  =  1  atm 
T  =  500®K 


50%  MASS  LOADING  OF 
PARTICLES  IN  JET 
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BUILDIIVIG-BLOCK  APPROACH  FOR 
ETC  GUN  SIMULATION 


1 .  NUMERICAL  STUDIES  WITH  EXISTING  VERSION  OF 
CRAFT 

PLASMA  JET-WORKING  FLUID  INTERACTIONS 
:  GRID  DYNAMICS/ADAPTIVE  GRIDDING 
:  THERMOCHEMICAL  EFFECTS 
:  BL/WALL  TEMPERATURES 
:  INJECTOR  EFFECTS/SINGLE  vs 
MULTIPLE 

;  STRONGLY-COUPLED  PROJECTILE 
SOLUTION 

2.  NEAR-TERM  UPGRADES 

-  COMPLEX  EQUATIONS  OF  STATE  -»•  Roe/TVD 
[RESEARCH] 

-  VOLUMETRIC  EFFECTS/DROPLETS 

-  INTERIOR  BALLISTIC/CAPILLARY  COUPLING 

-  DROPLET  FORMATION/COMBUSTION 
RADIATIVE  EFFECTS 

3.  PHENOMENOLOGICAL  STUDIES 

HOMOGENOUS  vs  NON-HOMOGENOUS 
MIXTURE  ASSUMPTIONS 

-  combustion  CHEMISTRY-SIMPLIFIED 
vs  MULTI-STEP  REPRESENTATION 

-  MIXING/PENETRATION  OF  PLASMA  JET- 
ke  vs  LES  TURBULENCE  MODELING 


SfetoiwAHpftMtionsiritBfTirtWfM/Cpypowtfoff 


EQUATION  OF  STATE  (EOS)  AND 
CHARACTERISTIC-BASED  UPWIND  SCHEMES 

•  RIEMANN  SOLVER  {Roe/TVD,  Godunov,  etc.) 
ORIGINALLY  FORMULATED  FOR  IDEAL  GAS  DYNAMICS 

•  MORE  RECENTLY,  ADVANCES  FOR  MULTI- 
COMPONENT,  REAL  GAS  FLOWS  WITH  CHEMICAL 
NONEQUILIBRIUM  (Vinokur,  Liu  -  1988-89) 

•  ETC  FLOWFIELDS  CHARACTERIZED  BY  VERY  COMPLEX 
EOS  WITH  COMPLICATED  CONSTITUTIVE  MODELING, 
PHASE  CHANGES,  etc.  (NON-CONVEX  EOS) 

-  EXACT  RIEMANN  SOLUTION  ENTAILS  PROHIBITIVE 
COSTS  FOR  LARGE-SCALE  COMPUTATIONS 

-  RECENTLY,  EFFICIENT  ALTERNATIVE  RIEMANN 
SOLVERS  FORMULATED  WHICH  COMPARE 
FAVORABLY  WITH  MORE  RIGOROUS  GODUNOV 
SCHEME 

GODUNOV  BASED  RIEMANN  SOLVER 


Collins,  P.  and  GIaz,  H.M.,  "Efficient  Solution  Algorithms  for  the 
Riemann  Problem  for  Real  Gases,"  J.  Computat.  Physics,  59  (1985), 
pp.  264-298. 

ALTERNATIVE  RIEMANN  SOLVERS 

Bell,  J.B.,  Coiella,  P.  and  Trangenstein,  J.A.,  "Higher-Order  Godunov 
Methods  for  General  Systems  of  Hyperbolic  Conservation  Laws, 

J.  Computat.  Physics,  82  (1989),  pp.  362-397. 

Davis,  S.F.,  "Simplified  Second-Order  Godunov  Type  Methods,"  SIAM 
J.  Sci.  Statit.  Comput.,  9  (1988)  pp.  445-473. 

Harten,  A.  Lax,  P.D.,  and  van-Leer,  8.,  "On  Upstream  Differencing  and 
Godunov-type  Schemes  for  Hyperbolic  Conservation  Laws,  SIAM, 

Rev.,  25  (1983)  pp.  35-61. 
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OBLIQUE  SELF-SIMILAR  SHOCK  REFLECTION 


F'lKuri*  1.  Si'lKuiiuUc  for  flowiield  luiliulizalion. 


NON-CONVEX  EOS 


CONVEX  EOS 


MODIFIED  POLYTROPIC 


POLYTROPIC 


p  =  \/-y  +  fiv) 
(PHASE  CHANGE) 


'•  An  Smokffee-Owneit  Company 


SeimKeAppiiGallamkmnmtionalCorpontion 
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p  =  V’y 


2D  Greenfarm  Simulation 


120miii 

Experimental  data  exhibits  pressure  waves 


VEL  (M/S) 


TEMP  (K)  „  PRESS  (PSI) 


Science  Applications  International  Corporation 
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TEMP  (K)  ^  PRESS  (PSI) 


Sdance  Applications  International  Corporation 
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An  Emptoyee-Onmed  Camper^ 


TEMP  (K)  „  PRESS  (PSI) 


Science  Applications  International  Corporation 

r  PRESSURE  VS  X  AT  0.39  MS 


DOOOO 


.0  p  temperature  vs  X  AT  0.39  MS 


10000.0 


1000.0 


100.0 


.2  0.2  ■ 

mm 

m- 
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.. ....  An  £m()io^9e-Ovime(i  Company 


TEMP  (K)  „  PRESS  (PSI) 


Science  Applications  International  Corporation 


PRESS  (PSI) 


BASE  PRESSURE  VS  TIME 


TIME  (MS) 


ePX-OFB 
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2D  Greenfarm  Simulation 

•  Model  provides  excellent  resolution 

•  Dynamic  gridding  &  moving  projectile 

•  Pressure  vi^aves  interact  with  geometry 

•  Complicated  flow  structure 


Summary 

CRAFT  applicable  to  ETC  flows 

-  state-of-  the-art  numerics 

-  synergistic  with  other  DoD  &  NASA  efforts 

-  emphasis  on  first-principles  model 

-  core  model  for  submodules  developed  by  DoD  & 
DOE  labs,  university,  contractors 

Preliminary  application  to  Greenfarm  data 

-  excellent  resolution  of  flow  structure 

-  geometric  effects  evident 

-  pressure  waves  observed 

Building-block  approach 

-  plasma-working  fluid  interaction 

-  liquid  EOS 


RECENT  ADVANCES  IN  CAP^^  GUN  MODELING 


Dave  Cook  and  Jahn  Dyvik 
FMC  Corporation 
Minneapolis,  MN  55421-1498 


ABSTRACT 


The  FMC  Corporation  is  continuing  the  effort  to  develop  and  adapt  ETC 
interior  ballistics  models  to  reflect  experimental  performance  of  present  gun 
designs.  A  new  version  of  a  2-D  Interior  ballistics  code  with  an  empirical 
combustion  model  has  been  applied  to  recently  obtained  30 -mm  data.  FMC  is  also 
testing  alternative  propellants  in  the  30-miii  gun,  and  efforts  to  describe  their 
thennochemistry  and  ballistic  performance  will  be  discussed. 

FMC  is  currently  firing  60-iBm  ETC  guns  under  contracts  to  the  Navy  and  DNA, 
as  well  as  under  IR&D.  These  efforts  are  testing  a  number  of  new  cartridge 
designs.  A  modeling  effort  has  recently  been  completed  which  describes  two  of 
these  new  designs  using  a  variable  area,  1-D  interior  ballistics  code.  This 
code  contains  the  same  gas  generation  model  as  the  2-D  code,  but  Includes 
additional  features  such  as  a  flexible  computational,  domain  and  an  empirical 
reaction  kinetics  model.  When  these  new  computational  algorithms  mature,  and 
cartridge  designs  are  decided,  the  new  1-D  features  will  be  incorporated  into 
the  2-D  code. 
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Recent  Advances  in  CAP.  Gun 

tm 

Modeling 


Dave  Cook 
Jahn  Dyvik 


Overview 


•  Status  of  2-D  Modeling  at  FMC 

•  1-D  Interior  Ballistics  Code  Development 

•  Recent  1-D  Applications  at  60-mm 

•  Impact  of  1-D  Calculations  on  Future  2-D 
Work 
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Status  of  2-D  Modeling  at  FMC 


The  2-D  model  of  the  CAP,^  process  contains  three  coupled  sub¬ 
models,  a  PFN  routine,  a  plasma  injector  routine,  and  a  2-D 
interior  ballistics  code  which  solves  the  viscous  Navier-Stokes 
equations. 

The  first  version  of  this  code  was  validated  using  mostly  large 
caliber  data.  Recentiy,  the  numerics  have  been  changed  to  a 
donor  cell  approach  and  an  artificial  diffusion  modei  has  been 
added. 

This  new  version  has  been  compared  to  30-mm  ballistics  data. 
Agreement  is  good,  however,  some  of  the  details  in  calcuiated 
pressure  traces  indicate  a  more  complex  combustion  model  is 
necessary. 


396 
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Empirical  Combustion  Model 


•  Three  functions  are  used  to  describe  the  energy  released  from  propellant 
combustion.  They  are: 

1)  the  gas  generation  ratethroughoutthechamber,  i.e.  the  "global"  rate  function, 

2)  the  spatial  distribution  of  the  gas  generation  rate,  and 

3)  An  Arrhenius  equation  for  the  gas  phase  reaction  rate 

•  The  "global"  rate  function  contains  two  parameters.  They  describe  the 
distribution  centroid  and  width.  These  are  the  most  critical  parameters 
and  must  be  chosen  by  considering  electrical  input,  propellant  properties, 
projectile  mass,  etc. 

•  The  spatial  distribution  function  is  calculated  by  considering  the 
normalized  product  of  two  densities.  The  product  function  of  gas  and 
liquid  density  is  likely  to  be  a  maximum  near  the  gas-liquid  interface.  It 
may  be  reasonable  to  assume  the  spatial  maximum  of  the  gas  generation 
rate  is  near  this  interface.  Two  exponential  parameters  are  used  to  shift 
the  product  maximum  if  desired. 
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DNA/SAIC  ETC  PROGRAM 


O 

o« 


30-inin  DNA  SHOT  2 


TIME  [mSl 


Thermochemistry 


•  2-D  30mm  and  1-D  60mm  calculations  have  been  performed 
using  peroxide/hydrocarbons  as  a  propellant. 

•  FMC  is  currently  testing  propellants  other  than  peroxide/ 
hydrocarbons  with  superior  properties  for  weapon  ization .  Tests 
include  propellants  consisting  of  HAN/hydrocarbons. 
Thermochemical  properties  have  been  determined  using  Blake. 

•  Thermochemical  studies  using  Blake  are  being  performed  to 
characterize  other  advanced  propellant  concepts  as  well. 
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Motivation  for  1-D  Code  Deveiopment 


•  FMC  employs  workstations  at  remote  testing  sites,  where  real  time 
execution  of  models  is  an  issue.  1-D  calculations  often  give 
adequate  representation  of  routine  ballistic  measurements  in  a 
fraction  of  2-D  execution  time. 

•  Relative  simplicity  of  1  -D  codes  allows  impact  of  cartridge  design 
details  on  ballistic  performance  to  be  estimated  rapidly.  A 
generalized  2-D  computational  domain  requires  much  greater 
effort. 

•  Theoretical  sub-models  can  be  tested  quickly  within  a  1-D 
environment.  Generalization  from  lower  to  higher  dimensionality 
is  often  easier  than  starting  from  scratch. 
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Formulation  of  1-D  Models 


•  The  l-D  models  are  collections  of  subroutines  that  describe  the 
PFN,  piasma  injector,  and  the  reacting  flow.  They  solve  1-D  Euler 
equations  with  artificiai  diffusion. 

•  A  donor  celi  method  is  used.  A  constant  number  of  ceils  of  equal 
thickness  expand  to  ailow  for  projectile  motion. 

•  The  cell  volume  and  face  area  are  design  dependent.  Cartridge 
features,  such  as  chamberage,  can  be  introduced  by  spatial 
dependence  of  cell  volume  and  face  area.  As  the  grid  expands, 
new  areas  and  volumes  are  calculated. 

•  Source  terms  in  the  Euler  equations  are  aiso  design  dependent. 
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1-D  Model  Setup  - 

Navy  60-mm  Cartridge  Design  A 


The  area  &  volume  of  each  computational  cell  has  been 
independently  calculated  to  reflect  the  following  design 
features: 

•  The  injector  flow  is  introduced  into  the  interior  of  the 
computational  domain;  flow  is  in  the  reverse  direction.  The 
empirical  propellant  -  plasma  model  has  been  changed  to 
account  for  this. 

•  The  domain  was  changed  to  describe  the  shape  of  the 
projectile  tail.  Algorithms  for  computation  of  the  base 
pressure  and  work  performed  were  also  changed.  Artificial 
diffusion  was  added  to  smooth  the  abrupt  change  in  area  at 
the  projectile  tail  interface. 


NJECTOR  PROPELLRHT  PROJECTILE 


PRESSURE  (  KPSI )  PRESSURE  (  KPSI ) 


0.0  1.0  2.0  3.0  4.0  5.0  6.0  7.0 

TIME  (  SEC  )  no"" 
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PRESSURE  (ksi) 

0.0  lao  20.0  30.0  40.0 
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1-D  Model  Setup  - 
DNA  60-mm  SFI  Design 


The  following  adaptations  have  been  made  to  the  baseline  algorithms: 

•  The  grid  extends  only  from  the  area  forward  of  the  fuel  rich  propellant 
to  the  rear  of  the  projectile.  A  regression  velocity  for  the  fuel  rich 
propellant  surface  is  calculated  so  that  the  injector  volume  is  added  to 
the  domain. 

•  The  resistance  of  the  injector  is  modeled  by  adjusting  the  geometry  of 
a  cylindrical  volume  to  match  experimental  results.  Mass  flow  from  a 
polyethylene  tube.  The  mass  flow  is  assumed  to  be  largely  fluid,  not 
gas. 
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DSHD 


SFI  -  ID  FLEXIBLE  GEOMETRY 
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60-mm  DNA  SHOT  #1 
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CAP  1-D  PRESSURE  DISTRIBUTION 

«  ^  L15  ihs  TRAVEL  =  7.06  cm 
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Impact  of  1-D  Calculations  on  2-D  Work 


•  Effects  due  to  cartridge  geometry,  i.e.  chamberage,  projectile 
tail  boom,  etc.,  were  found  to  have  significant  impact  on 
ballistic  calculations.  The  2-D  computational  domain  must 
be  generalized  to  better  reflect  the  hardware. 

•  The  species  equations  and  reaction  kinetics  model  first 
used  in  the  60-mm  DNA  SFI  calculations  were  found  to  have 
several  beneficial  effects.  For  instance,  spatial  variation  in 
chamber  o/f  cannot  be  included  without  tracking  species 
equations.  Also,  the  "lag"  time  that  the  Arrhenius  equation 
introduces  between  the  chemical  energy  release  rate  and 
gas  generation  rate  inhibited  the  fuel-oxidizer  reaction  during 
electrical  pulse  onset. 
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aO-MM  ETC  BALLISTIC  DIAGNOSTIC  FACILITY 


K.  White,  I.  Stobie,  B.  Benslnger,  S.  Driesen, 
H.  Burden,  G.  Katulka,  and  A.  Juhasz 
US  Army  Ballistic  Research  Laboratory 
Aberdeen  Proving  Ground,  HD  21005-5066 


ABSTRACT 

The  BRL  has  installed  a  BO-nrni  Electrothermal  Chemical  ballistic  diagnostic 
facility,  built  under  contract  by  the  FMC  Corporation.  The  fixture  has  the 
capability  of  goin^  to  pressures  as  high  as  650  HPa  with  a  maximum  chamber  volume 
of  550  cc  and  a  projectile  travel  of  1.53  m  (50  caliber) .  Pressure  and  temperature 
measurements  can  be  made  in  the  chamber  and  down  bore.  The  fixture  was  built 
in  modular  form  such  that  the  conventional  steel  chamber  can  be  replaced  by 
either  an  optically  or  x-ray  transparent  chamber,  allowing  for  flow  visualization 
during  the  electrothermal  process.  In  addition  to  the  gun  mode,  the  fixture  can 
be  operated  in  either  the  closed  chamber  or  burst  disc  mode  to  examine  the  effect 
of  a  static  boundary  condition  on  the  ETC  process. 

The  fixture  was  test  fired  three  times  in  the  gun  mode  with  projectile 
velocities  of  over  1000  m/s,  peak  pressures  of  280  HPa  and  with  60  kJ  of  electrical 
energy  input.  Two  tests  were  performed  tin  the  burst  disc  mode  (65  HPa). 

Tests  will  also  be  conducted  with  water  as  a  working  fluid  and  with  various 
plasma-water  Interface  geometries.  These  tests  are  to  act  as  a  baseline  test 
series  of  plasma-fluid  Interactions  without  chemical  contribution  from  the  working 
fluids.  Pressure  waves  have  been  observed  in  ETC  firings  with  reactive  fluids 
and  these  tests  will  help  to  isolate  the  hydrodynamic  processes  from  the  combustion 

?rocesses.  These  tests  will  be  carried  out  with  projectiles  (moving  Boundary 
ayers)  and  blow-out  discs  to  Isolate  the  effect  of  the  moving  boun^ry  layer 
on  the  mixing  process.  In-bore  projectile  displacement  measurements  will  also 
be  carried  out.  Tests  are  in  progress  and  will  be  discussed. 
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LABORATORY  COMMAND 

BALUSTIC  RESEARCH  LABORATORY 


30-MM  ETC  BALLISTIC 
DIAGNOSTIC  FACILITY 


K» 

K» 


K.  White,  I.  Stobie,  B.  Bensinger,  S.  Driesen 
H .  Burden ,  G .  Katulka  and  A .  Juhasz 


JANNAF  Workshop 
ETC  Modeling  &  Diagnostics 
July  9-11,  1991 
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US  ARMY''^^ 
LABORATORY  COMMAND 


PROGRAM  OBJECTIVE 


Develop  data  base  of  ETC  firings  for  model  validation: 

-  Axial  and  radial  plasma  injection  configurations , 

-  Inert  and  reactive  working  fluids , 

-  Moving  boundary  layer  (projectile)  &  static , 

-  Short  &  long  electrical  pulse , 

-  Optically /x-ray  transparent  chamber. 
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BALLISTIC  RESEARCH  LABORATORY 


30-mm  fixture 

maximum  chamber  volume ,  550  cc , 
maximum  operating  pressure ,  650  MPa , 
projectile  travel,  1 .53  m  (50  caliber), 
gun  or  blow-out  disc  mode , 
measurements , 

pressure ,  chamber  and  barrel , 

projectile  displacement ,  35  GHz  interferometer , 
resolution  1mm, 

optical/x-eray  access  chamber. 
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TEST  CONDITIONS 


US  ARMT 
LABORATORY  COMMAND 


BALLISTIC  RESEARCH  LABORATORY 


Scaling:  (30/120)®  =  1/64 


120-mm 

30-mm 

Electrical  energy 

4 . 5  MJ 

70  kJ 

Chamber  volume 

9.95  1 

155  CC 

Projectile  mass 

11.4  kg 

178  g 

Problem  with  scaling  for  hydrodynamic  effects . 


Axial  and  radial  plasma  injection: 

-  Inert  working  fluid,  water, 

-  Moving  boundary  layer ,  projectile , 

-  Short  electrical  pulse ,  300  [is ,  70  k J . 
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INJECTOR  & 


BALLISTIC  RESEARCH  LABORATORY 


Radta)  (Center  Core) 


us  ARMY 

LABORATORY  COMMAND 
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Radial  (  Center  Core  ) 


Modified  BRL  five  module  pulse  forming  network  (135kj  at  10k V) 
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PRESSURE  GAGE  LOCATION 


BALLISTIC  RESEARCH  LABORATORY 


US  ARMY 
LABORATORY  COMMAND 


Rll  dimensions  in  mm  . 


3 

Chamber  Uolume  »  547.4  cm 
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FOURIER  ANALYSIS 
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BALLISTIC  RESEARCH  LABORATORY 
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FOURIER  ANALYSIS 


US  ARMY  >5^ 
LABORATORY  COMMAND 


BALLISTIC  RESEARCH  LABORATORY 

Radial  Axial 

W/O  Liner  W  Liner 

070291  061991 


Axial 
W/O  Liner 
062091 


BALU8TIC  RESEARCH  LABORATORY 


Radial  Injection  (Center  Core) 


ID 

P2 
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V 
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EE 
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kJ 
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Note:  Water  leakage  around  projectile  in  test  ID  062591 . 
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CONCLUSIONS 

LABORATORY  COMMAND 

BALLISTIC  RESEARCH  LABORATORY 

Radial  (Center  Core)  vs .  Axial  Injection 
Radial 

-  more  gradual  pressure  rise  at  projectile , 

-  lower  pressures  and  velocities , 

-  pressure  drops  off  more  rapidly . 

Issues: 

-  The  water/plasma  energy  system  does  not  behave 

as  a  well  stirred  system . 

-  Pi  levels  suspect , 

-  Chamber  liner  may  attenuate  high  frequency  waves 

acoustic  attenuation, 

longer  distance  from  pressure  to  gage , 

(10-15  mm). 

-  Pressures  contain  substantial  high  frequency  components , 

not  combustion  driven , 

PROPELLANT  SURFACE  AREA  IMPLICATIONS . 


BALUSnC  REBEABCH  tiABORATORT 


Deyelop  data  base  of  ETC  firings  for  model  validation; 

-  Axial  and  radial  plasma  injection  configurations , 

-  Inert  and  reactive  working  fluids , 

-  Moving  boundary  layer  (projectile)  &  static , 

-  Short  &  long  electrical  pulse 

-  Optically /x-ray  transparent  chamber , 

Reactive  working  fluid: 

la-molar  HAN .  NOS-365 ,  HAN  +  IP  AN  +H20 , 

low  energy ,  safe  &  available , 

safe  &  available ,  difficult  to  ignite , 

some  ballistic  data .  bulk  loaded  firing  data , 

mono-propellant . 


VOLTS 


INTERFEROMETER  DATA 
&  DISPLACEMENT 


utl  efal991.  l«TEKF£ROnETER 


TINE 
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NUMERICAL  SIMULATION  OF  THE  INTERIOR  BALLISTIC  PROCESSES 

IN  AN  ETC  GUN 

J.  L.  Chen,  F.  B.  Cheung  and  K.  K.  Kuo 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


ABSTRACT 

A  comprehensive  theoretical  model  has  been  developed  to  study  the  Interior 
ballistic  processes  in  an  ETC  gun.  Phenomena  occurring  in  five  separate  regions 
are  considered  In  the  model.  These  regions  consist  of  the  gas  phase  (l.e., 
Taylor  cavity),  working  fluid  (l.e.,  liquid  propellant),  dispersed  droplets  in 
the  Taylor  cavity,  projectile,  and  plasma  generating  cartridge  (PGC).  Dispersed 
droplets  are  generated  through  the  entrainment  process  as  a  result  of  the  relative 
motions  between  the  gas  and  llauld  phases.  Governing  equations  are  formulated 
from  the  first  principles  of  physics,  accounting  for  the  plasma/working  fluid 
Interactions  and  the  projectile  dynamics.  TWo  different  sets  of  PFN  discharge 
curves  have  been  employed  to  identify  the  effects  of  PFN  discharge  characteristics 
on  gun  performance.  Various  energy  ratios  (chemical  energy /electrical  energy) 
have  been  used  to  Investigate  how  the  propulsive  ener^  is  Influenced  by  plasma 
generation  and  combustion  of  liquid  propellant  in  different  stages  of  the  ETC 
gun  event.  The  model  has  been  partially  validated  by  comparing  the  computed 
breech  pressure-time  traces  with  the  measured  data  provided  by  FHC. 
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OBJECTIVE 


*  Develop  a  Comprehensive  Theoreticar  Model  to  Investigate 

the  Interior  Ballistic  Processes  of  an  Electrothermal 
(ET)  Gun 

*  Predict  the  Interior  Ballistic  Processes  of  an  ET  Gun 

*  Seek  a  Better  Understanding  of  the  Mechanisms 

Controlling  the  Performance  of  an  ET  gun 
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Working  Fluid 


d.  Intensive  Mixing  and  Combustion 


Fig.  1-2  Schematic  of  a  Typical  ET  Gun  Firing  Cycle 
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Control  Volume  for  Gas  Motion  of  Droplet 

=  Total  Control  Volume 
-Volume  Occupied  by  Droplets 


Control  Volume  for 
Liquid  Phase 


Projectile  Dynamics 


Control  Voluntics  for  Different  Phases  in  an  ET  Gun 


THEORETICAL  FORMULATION 
Assumptions 

1.  Gas.  Liquid,  and  Solid  Phases: 

Transient  Axisymmetric  Cylindrical  Eulerian 
Coordinate 

2.  Motions  of  the  Droplets  : 

Small  Droplets: 

*  Following  Gas  Flow  Motion 

*  Number  Density  Equation 

3.  Combustion  of  Liquid  Propellant  : 

Oxidizer  Droplets  Burning  in  the  Gaseous  Fuel 
Environment 

(Heat  of  Vaporization  of  the  Oxidizer  is  Much  Higher 
than  that  of  the  Fuel  at  Standard  State) 

4.  Combustion  Mechanisms  of  Liquid  Propellant  : 
Simple  Chemical  Reacting  System  (SCRS) 

*  Oxidizer,  Fuel,  and  Product 

*  One  Step  Chemical  Reaction 

5.  Body  Force  Effect  :  Negligible 
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Governing  Equations  for  Gas  Phase 

*  Deriving  from  First  Principles  of  Physics 

*  Continuity,  x-,  r-Momentums,  and  Energy  Equations 


*  Turbulence  Closure 

1.  Favre  Averaged 

2.  K,  £  Two  Equation  Model 
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Gas  Phase 


Continuity  Equation 
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Gas  Phase 


Energy  Equation 
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Equation  of  State 


Gas  Phase 


*  Turbulence  Model 
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Droplet  Phase 


Droplet  Motion: 

Small  Droplets  Follow  Flow  Motions 
Number  Density  Equation  for  Group  k 


d[n^  a(nkUg)  aCn^Vgr) 

at  ^  a"x  ^  ra r 


M*1 


(a)  (b)  (c) 


(a) :  Due  to  Combustion  of  Group  k 

(b) :  Due  to  Combustion  of  Larger  Droplets 

(c) :  Due  to  Entrainment  Mechanisms 


Droplet  Combustion  : 

Boundary  Condition  at  Droplet  Surface 
Pi,g  =  Pi,l  (at  Ambient  Pressure) 
fi,g  =  fid  (at  High  Pressure),  Equality  of  Fugacity 
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Entrainment  Mechanism 


/ 

Roll  Wave,  Wave  Undercut,  Bubble  Burst,  Liquid 
Impingement,  Liquid  Bulge  Disintegration 


SAS 


'i  }  .  i  )  )  }  )  )  i  , ,  ;  i 


TYPE  1 
ROLL  WAVE 


rrPE  3 

BUBBLE  BURST 


TTTTTrTTTTTTTTTTTTTTrrTTTTT 


TYPE  4 

UQUIO  IMPINGEMENT 


TYPE  5 
UQUIO  BULGE 
DISINTEGRATION 
(COUNTER  CURRENT) 
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Droplet  Phase 


Entrainment  Rate: 


Ra  =  K. 


Wlf-Wlf.c 


pf  pr^ 


Ra:  entrainment  rate 

Ka:  entrainment  coefficient 

Pe:  wetted  perimeter 

Wlf:  liquid  film  flow  rate 

Wlf.c:  critical  liquid  film  flow  rate 


Entrainment  Droplet  Size 

Average  Droplet  Size  Based  on  a  Log  Normal 
Distribution  for  Annular  Flow 


d32  f^S  *^8  dt 
di  I  2  a 


1/2 

=  2.3*  10 


d32:  Sauter  mean  diameter 
dt:  hydraulic  diameter  of  tube 
fs:  friction  factor  for  smooth  pipe 
a:  surface  tension 


155 


Projectile  Dynamic  Motions 


Projectile  Dynamic  Motions  : 

Sliding,  Rotating,  and  Engraving  in  the  Gun  Barrel 

Force  Balance  : 

Driving  Force  Pressure  Acting  on  Projectile  Base 
Resistant  Forces  :  Due  to  Sliding,  Rotating,  and 
Engraving  Motions 


.-^Numerical  Solutions 


*  Computational  Domain 


Gas  Cavity 


(c)  Annual  Liquid  Film  Stage 


(d)  Free  Expansion  stage 
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Coordinate  Transformation 


Original  Coordinate  :  t,  x,  r 


New  Coordinate  :  x, 


r  =  t.  T)=x4-Uc(t-to),  C  =  r  +  v^(t-to) 


a  (<i>vpg)  a  (ovpg) . ..  5(<i>Wg) . ..  ^  (‘i>¥pg) 

3t  3t  8ti 
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Algorithm 


*  SIMPLE  (Semi-Implicit  Method  for  Pressure-Linked 
Equations) 

*  Highly  Compressible  Flows  : 

a.  Density  Variation  Effects  on  Pressure 
Correction  Equations 

b.  Equation  of  State  -  Noble  Abel  Equation 


^Pressure  Correction  Equation 
^pPp=  aEPE+  a^Pw'*'  as^s 


aE=(Pe^e-KUe)ATi 
aw=(Pwdw-Kuw)  Ati 


as=(p*d3-Ku;)A; 

aN=(Pndn-Ku;)A^ 

ap=(P*d„+p‘d3)Aj;+(p;d,+  p^d^)  Ati  + 


A: 

dj=^'  J  =  e,w.s,n 
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RT 
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Peak  Power  Time/ 

Total  Discharging  Time 


PRESSURE  (PA) 


Breech  and  Taylor  Cavity  Front  Pressure-Time  Traces  with  PFN 
Discharging  Time  of  2.0  ms  and  Peak  Power  Time  at  0.8  ms 
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Pressure  (PA) 

4.3  E+8 

3.3  E+8 

2.3  E+8 

1.3  E+8 

0.3  E+8 


Time  History  of  Averaged  Axial  Pressure  in  the  ET  Gun  Chamber  with  PEN 
Discharging  Time  of  2.0  ms  and  Peak  Power  Time  at  0.8  ms 
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Pressure  (PA) 
2.78  E+8 


2.73  E+8 


2.68  E+8 


2,63  E+8 


2.58  E+8 


Spatial  Pressure  Distribution  in  the  ET  Gun  Chamber  at  0.33  ms  with  PEN 
Discharging  Time  of  2.0  ms  and  Peak  Power  Time  at  0.8  ms 


VELOCITY  (m/S) 


Projectile  and  Taylor  Cavity  Frontal  Velocity-Time  Traces  with  PFN 
Discharging  Time  of  2.0  ms  and  Peak  Power  Time  at  0.8  ms 
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o^ 

Oi 


Velocity  (pi/s) 
600  A 


0.013 


008 


r(m) 


x(m) 


Spatial  Distribution  of  Axial  Gas  Velocity  in  the  ET  Gun  Chamber  at  0.45  ms 
with  PFN  Discharging  Time  of  2.0  ms  and  Peak.Power  Time  at  0.8  ms 
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0.06 


0.013 


0.000 


0.00 


x(m) 


Spatial  Distribution  of  Radial  Gas  Velocity  in  the  ET  Gun  Chamber  at  0.45  ms 
with  PFN  Discharging  Time  of  2.0  ms  and  Peak  Power  Time  at  0.8  ms 


PROJ  VEL  (m/S) 


PROJ  VEL 
PROJ-TRAJ 


Projectile  Trajectory  and  Velocity-Time  History  with  PFN  Discharging  Time 
of  2.0  ms  and  Peak  Power  Time  at  0.8  ms 
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PROJ-TRAJ  (m) 


.5 


2 


1.0 


1.0  1.5  2.0 

Tfme  (ms) 


Breech  and  Taylor  Cavity  Front  Pressure-Time  Traces  with  PFN  Peak 
Power  Time  at  1.4  ms  and  Discharging  Time  of  1.6  ms 


POWER  (MW)  POWER  (MW) 


(a)  Various  Peak  Power  Times 


TIME  (mS) 


(b)  Various  Total  Discharging  Times 


TIME  (mS) 


PFN  Discharging  Power  Curves  with  (a)  Various  Peak  Power  Times  and 
(b)  Various  Total  Discharging  Times 


Results  from  Numerical  Analysis  for  Cases  with 
Various  Peak  Power  Times 


Interior  Ballistic  Data 

Peak  Power 

Time  (ms) 

.  Total 
Discharging 
Time  (ms) 

Maximum 
Breech 
Pressure  (PA) 

Muzzle 

Velocity  (m/s) 

0.8 

2.0 

2.11 

4.10*108 

1936 

12.0  % 

1.0 

'2.0 

■  2.18 

3.95*108 

'1924 

11.8  % 

1.4 

2.0 

2.32 

3.40*108 

1882 

11.0  % 

1.8 

2.0 

2.43 

3.39*108 

'1822 

10.6  % 

♦Electrical  Energy  0.7  MJ,  Chemical  Energy  0.4  MJ 


Results  from  Numerical  Analysis  for  Cases  with 
■'  ■  Total  Discharging  Times  "  ' 


PFN  Discharging 
Characteristics 

Interior  Ballistic  Data 

Peak  Power 

Time  (ms) 

Total 

Discharging 
Time  (ms) 

Maximum 
Breech 
Pressure  (PA) 

Muzirie 

Velocity  (m/s) 

Ballistic 

Efficiency 

1.2 

1.4 

2.08 

3.60*108 

2048 

13.3  % 

1.4 

1.6 

2.22 

3.40*108 

1972 

12.4  % 

1.8 

2.0 

2.43 

3.39*108 

1822 

10.6  % 

♦Electrical  Energy  0.7  MJ,  Chemical  Energy  0.4  MJ 


Interior  Ballistic  Results  from  Present  Work  and  Literature 


(a)  Measured  Results 


Gun  (mm) 

Barrel 

(caliber) 

Projectile 
Mass  (gm) 

Muzzle  Vel 
<km/S) 

Electrical  Eng 
(kl) 

Working 

Fluid 

BSSSI 

Max  Breech  P 
(psi) 

Greig  et  al 
(19881 

20 

mm 

8 

3.5 

* 

1% 

♦ 

* 

Greig  et  al 
(1988) 

20 

21 

2.8 

700 

« 

11.8% 

laSfM 

Greig  et  al 
(1988) 

20 

50 

2.1 

HI 

* 

15.7% 

* 

* 

Greig  et  al 
(1988) 

mm 

102 

1.43 

« 

5.2% 

*  ■ 

*  • 

10 

* 

1-6.5 

1.2  -  1.7 

1^-325 

exothermic 

propellant 

* 

Chryssom* 
alUs  (1988) 

■■ 

wm 

50  -  380 

1.1  -  2.75 

exothermic 

propellant 

mm 

* 

Chryssom* 
allis  (1988) 

mnn 

* 

mm 

hb 

exothermic 

propellant 

* 

Chryssom- 
allis  (1988) 

■■ 

* 

2000- 

6000 

0.7  -  1.5 

mui 

* 

■rJiMi'ilM 

rrr*rW 

*  Data  Unavailable  -fElectrical  Efficiency  »  Projeccile  Kinetic  Energy  /  Input  Elecuica!  Energy 


Ballistic  Efficiency  *  Projectile  Kinetic  Energy /Total  Input  Energy 


(b)  Predicted  Results 


Source 

Gun  (mm) 

Barrel 

(caliber) 

Projectile 
Mass  (sm) 

Muzzle  Vel 
(km/s) 

Working 

Fluid 

rnESmSm 

WKSam 

rnsmBsmi 

mnW 

Oberle  (1988) 

14 

18 

2.17 

H2O 

447.8 

Bpfff 

9% 

mhvim 

Obcrle  (1988) 

14 

100 

18 

2.34 

LiBHd 

449.8 

406.4 

■hmm 

mmm 

Oberle  (1988) 

14 

18 

1.75 

TiH2/Al 

77.9 

HMTiM 

12.5% 

Oberle 

(1988) 

14 

mm 

18 

2.13 

CgHis/ 

H7O2 

43.3 

322.6 

HI 

mmn 

Present 

Work-" 

25 

■■ 

70 

1.94  -  2.05 

CrHis/ 

H902 

700 

mmn 

mirnkm 

+  Max  Breech  Pressure  =  52,900  -  60,200  psi,  Cun  System  Information:  Projectile  Mass  =  70  gni.  Residual  Air  =  5%,  Chamber  Lengih/DLameicr  =  3.5. 

Electrical  Energy  =  0.7  MJ.  Total  Energy  =  1. 1  MJ,  Working  Fluid  =  C8H}8/H2C>2 


Conclusions 


1.  A  comprehensive  theoretical  model  describing  the 
interior  ballistic  processes  of  an  ET  gun:  has  been 
formulated. 

2.  The  effects  of  PEN  discharging  characteristics  on  gun 
performance  have  been  identified. 

3.  Muzzle  velocity  can  be  optimized  by  tailoring  the  PEN 
discharging  time  and  peak  power  time 

while  holding  a  moderate  maximum  breech  pressure. 
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INTENTIONALLY  LEFT  BLANK. 
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FINITE-ELEMENT  MODELING  OF  ELECTROTHERMAL-CHEMICAL  GUNS 


N.  K.  Wlnsor  and  S.  A.  Goldstein 
GT-Devices,  inc. 
Alexandria,  VA  22312 


ABSTRACT 

Two-dimensional  modeling  of  reactive  flow  dynamic 
cartridge  Is  described  and  Illustrated.  Special  emphasl 
effects  of  ullage  and  rate -dependent  chemistry  on  large- 


s  in  an  electrothermal 
s  will  be  placed  on  the 
-scale  wave  dynamics. 


177 


1991  JANNAF  ETC  Modeling  Workshop 

FINITE-ELEMENT  MODELING  OF 
ELECTROTHERMAL-CHEMICAL  GUNS 

N.  K.  WINSOR  &  S.  A.  GOLDSTEIN 
GT-DEVICES,  INC. 

A  subsidiary  of  Generai  Dynamics  Land  Systems 

57 05 A  General  Washington  Drive 
Alexandria,  VA  22312-2408 

703-642-8150 


SPECIAL  THANKS  TO: 

Simon  Wang 
Bob  Greig 
Hugh  McElroy 
Steven  Bunte 
William  Oberle 


stub  Case  Area 


C5ENERAL  DYNAMICS 
Land  Systems  Division 


GT  -  Devices 

CARTRIDGE  GEOMETRY 


Cartridge  Case 


Obturator 
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GT  -  Devices 


GENERAL  DYNAMICS 

Land  Systems  Division 


OUTLINE 

FINITE-ELEMENT  MODEL 
WATER  HAMMER  TESTS 
CLOSED  BOMB  EXPTS. 
IBHVG2  #6  SIMULATION 
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•  C3ENERAL  DYNAMICS 

eviC6S  Land  Systems  Division 

SUMMARY 

1 .  WATER  HAMMER:  SOUND  SPEED 
DEPENDENCE  ON  PRESSURE  IS  KEY 
TO  DYNAMICS  OF  ELECTROTHERMAL 
AND  LIQUID  PROPELLANT  GUNS. 

2.  CLOSED  BOMB:  AN  IDEAL  METHOD 
FOR  VALIDATING  BURN  MODELS. 

3.  IBHVG2  TEST:  UNSTRUCTURED 
FINITE-ELEMENT  METHODS  WORK 
SPECTACULARLY  IN  GUN  GEOMETRY! 


GT  -  Devices 


GENERAL  DYNAMICS 
Land  System's  Division 


FINITE-ELEMENT 
ETC  MODEL  (FETC) 

NAVIER-STOKES  EQ.  OF  MOTION 
4th  TIME,  2nd  SPACE  FCT  ALG. 
UNSTRUCTURED  FINITE-ELT.  GRID 
DETAILED  CARTRIDGE  GEOMETRY 
DETAILED  PROPELLANT  CHEMISTRY 


GT  -  Devices 


GENERAL  DYNAMICS 

Land  Systems  Division 


NAVIER-STOKES 


TRANSPORT  EOUATIONS 


alj/ :  dp;  i  dp: 

- 7“  “f"  ^ - ■■  ^ - 
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STATE  VECTOR  AND  FtUXES 
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GT 


Devices 


GENERAL  DYNAMICS 

Land  SY'stems  Division 


SCHEMATIC  GEOMETRY 


GT  -  Devices 


general  dynamics 

Land  Systems  Division 


BASIC  GRID  STRUCTURE 


GENERAL  DYNAMICS 
Land  Systems  Division 


SELECTIVELY-REFINED  GRID 
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GT  -  Devices 


GENERAL  DYNAMICS 

Land  Systems  Division 


WATER  HAMMER  TESTS 

WATER  EQUATION  OF  STATE 
COMBUSTIBLE  STUB  CASE  "VOID" 
PRESSURE  WAVE  ->  SHOCK  WAVE 
LARGE  AMPLITUDE  TRANSIENTS 
HOW  TO  CONTROL  THEM 


Sound  Speed  (km/sec) 


GT 


Devices 


G 


SOUND  SPEED  VS.  TEMPERAI 
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GT  -  Devices 


GENERAL  DYNAMICS 

Land  Systems  Division 


120mm  GEOMETRY 


COMBUSTIBLE  CASE,  VIEWED  AS 


INITIAL  VOID,  LIMITING  SOUND  SPEED. 
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GT  -  Devices  cseneral  dynamics 

— -  Land  Systems  Division 

SOUND  SPEED 

REDUCED  BY: 

LOW  PRESSURE 
BUBBLES 
OTHER  VOIDS 

EFFECTS: 

SNOWPLOW  SHOCK 
WATER  HAMMER 


GT-Devices 


GENERAL  DYNAMICS 

LMd  SyBtema  Division 


CASE  1.  UNIFORM  POWER  INPUT 


POWER  VS.  TIME 


5 


H- 

3- 

2- 


1 


0  too  200  300  400 


Time,  microseconds 


PRESSURE  VS.  TIME 


1  -  IGNITER  2  -  REFLECTOR 


3eA602^ 


Power.  MW/cc 


GT-Oeviees 

CASE  2.  TAILORED 

POWER  VS.  TIME 

5 

4 

3 

2 

1 

0 


0  too  200  300  400 

Time,  microseconds 


re.  KRobore 
Ajeonds) 


J6A6022 
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CaENERAL  DYNAMICS 
Land  Systems  Division 


GT  -  Devices 

CONTROL  OF  TRANSIENTS 

SEND  A  COMPRESSION  WAVE 
MINIMIZE  THE  SNOWPLOW 
MEET  THE  REFLECTION 


BE  GENTLE! 


GENERAL  DYNAMICS 
Land  Systems  Division- 


GT  -  Devices 

CLOSED-BOMB  EXPTS. 

FIXED-WALL  GEOMETRY 
FULL  MODEL  CHEMISTRY 
TRANSIENT  COMBUSTION 


MOVING-WALL  GENERALIZATION 


GT  -  Devices 


general  dynamics 

Land  Systems  Division 


VARIOUS  CLOSED-BOMB  CURVES 

dP/dt  VERSUS  P 


C'osed— bomb  curve  Ctos«d-'bomb  curve  Closed  -  t>c>n>b  curve  Closed-bomb  curve 


Closed— bomb  curve  Ciosed-bomb  curve  Closed-*bumb  curve  Closed— bori>b  curve 


«tp« 
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GT  -  Devices 


GENERAL  DYNAMICS 

Land  Systems  Division 


IBHVG2  TEST  CASE  #6 


FEATURES  INCLUDED: 

M829  LONG  ROD  PROJECTILE 
Ml  25  PRIMER  TUBE  FLAME  SPREAD 
JA2  COMBUSTION  RATE 

FEATURES  OMITTED: 

PROPELLANT  PERF  GEOMETRY 
JA2  COMPRESSIBILITY  IN  EOS 

RESULT: 

RELATIVELY  COMPLETE  SIMULATION 


RftTE  ( IN/S) 


Temperature  Sensitivity 

US  ARMY 


LABORATORY  COMMAND 


BALLISTIC  RJESEARCH  LABORATORY 


Solid  Propellant  JA2 
Experimental  closed  chamber  burn  rates 
at  Ambient,  at  5Q°C,  and«32°C 


DISTANCE  (M) 


GT 


general  dynamics 

Land  Systems  Division 

1  20  MM  SIMULATIONS  COMPARED 

_  IBHVG2  AND  FETC  PROJECTILE  TRAVEL 

4 

3.5 

3 

2.5 

o 

Am 

1.5 

1 

0.5 

0 


0  2  4  5 

TIME  CMS) 


Devices 


PRESSURE  (KBAR) 


GT  -  Devices 


GEMERAL  DYNAMICS 

Land  Systems  Division 


120  MM  IBHVG2  SIMULATION 


STUB  AND  PROJECTILE  BASE  PRESSURES 


TIME  (MS) 


-Devices 


CaENERAL.  DYNAMICS 

Land  Systems  Division 


120  MM  FETC  SiMULATION 

STUB  AND  '  PROJECTILE  BASE  PRESSURES 


0 


TIrrw,  mllllW'”<3rKte 


Pressure  (Atmospheres) 
(Thouson<ls) 


GT 


CaENERAL  DYNAMICS 

Land  Systems  Division 

120  MM  FETC  SIMULATION 

_  DETAIL  OF  PROJECTILE  BASE  PRESSURE 

4 

3.  a 

.5 

2.6 

2 

1.5 

1 

0.5 

0 

Tims. 


Devices 
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GENERAL  DYNAMICS 

Land  Systems  Division 


GT  -  Devices 


120  MM  FETC 
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general  dynamics 

Lsnd  Systems  Division 

LAYOUT 


_  _■  .  cbcip  GENERAL  DYNAMICS 

GT  -  pGVICeS  .  ..  Systems  Division 

SUMMARY 


1 .  WATER  HAMMER:  SOUND  SPEED 
DEPENDENCE  ON  PRESSURE  IS  KEY 
TO  DYNAMICS  OF  ELECTROTHERMAL 
AND  LIQUID  PROPELLANT  GUNS. 

2.  CLOSED  BOMB:  AN  IDEAL  METHOD 
FOR  VALIDATING  BURN  MODELS. 

3.  IBHVG2  TEST:  UNSTRUCTURED 
FINITE-ELEMENT  METHODS  WORK 
SPECTACULARLY  IN  GUN  GEOMETRY! 


Intentionally  left  blank. 
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SPECIAL  DIAGNOSTICS  AND  INSTRUMENTATION 


Rex  D.  Richardson  and  Michael  D.  Haworth 
Science  Applications  International  Corporation 
Albuquerque,  NM  87106 


ABSTRACT 

SAIC  and  FMC  experimenters  are  currently  planning  to  field  two  special 
diagnostics  for  application  on  the  Defense  Nuclear  Agency's  ETC  Propulsion  for 
Enhanced  Gun  System  Effectiveness  program.  A  three-frame  flash  x-ray  radiographic 
system  Is  being  developed  to  observe  dynamic  mixing  processes  In  the  30  mm  FMC 
CAP^  test  fixture.  Results  of  x-ray  transmission  calculations  along  with  cold 
mod^  flash  x-ray  tests  at  300  kV  will  be  presented.  Data  analysis  and  Image 
processing  techniques  will  be  described.  The  second  diagnostic  Is  a  laser 
rangefinder  system  designed  to  track  the  position  of  a  projectile  In-bore  from 
shot  start  to  muzzle  exit  at  a  resolution  of  less  than  one  centimeter  and  a  data 
rate  of  AO  kHz.  System  design  Issues  and  the  prototype  test  plan  will  be 
discussed. 
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JANNAF  Workshop  on 
ETC  Modeling  and  Diagnostics 

June  1991 


dlfftlCD 

Defense  Nuclear  Agency 


SPECIAL  DIAGNOSTICS  AND 
INSTRUMENTATION"' 


Prepared  by 


Dr.  Rex  D.  Richardson 


Science  Applications  International  Corporation 
Albuquerque,  New  Mexico 

cooiVte.  sui>noit.rx:i>  S-y  Ths 


•PMC 


'  Science  Appileatlone 
International  Corpt^adon 

ArEirphjm-OmmlCemrMif 
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INTRODUCTION 


•  X-Ray  radiography  of  fluid/plasma  jets 

Overview 

Technical  approach 
FSCATT  runs 
300  kV  FXR  tests 
Image  analysis 

•  Laser  rangefinder  in-bore  diagnostic 

Conceptual  design, system  issues 


Sc/ence  AppUcatlona 
Int9matton0l  Corporation 

An  £iivlo)mn-O$n0ii  Conpanf 


Rn  81/680  ha/da 


FLASH  X-RAY  OBSERVABLES 
CAP  SYSTEMS 


COMBUSTION  CHAMBER 


TAYLOR  CAVITY 

•  Volume 

•  Density 


TURBULENT  ZONE 
ENTRAINMENT  &  MIXING 

JET  FRONT 
PROPAGATION 


Science  Applications 
International  Corporation 

An  fir^|wOim0tf  CofT^/iy 


RR  91/797  rvs>da 


p  Sdenea  Applications 
7^  tntamation^  Corporation 

AnEmplejVt  ftwxtfCwnpw 


RR  91/860  hWda 


LAYOUT  OF  3  FRAME  FXR  SYSTEM 
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FSCATT  CODE  RUNS 


•  X-ray  transport  and  energy  deposition 

•  Discrete  ordinate,  ID 

•  PE  absorption,  Compton  scattering  and  fluorescence 

•  Input/output  energy  flux  and  fliience 

•  Angular  distribution  of  secondary  radiation 


X-RAY  TRANSMISSION  MODELING 
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FSCATT  INPUT  GEOMETRIES 


Abbreviations 


cc 

Carbon  composite 

w 

Tungsten 

ox 

Oxidizer 

Ba 

Barium 

F 

Fuel 

Au 

Gold 

Mp 

Mass  fraction  of  High  Z 
dopant  in  Fuel 

WIqx 

Mass  fraction  of  High 
Z  dopant  in  Oxidizer 

Case  1 


layer 

thickne$s 

(inches) 


- > 

2t 


1  =  0.5,0:75,1.0 

-f  ' 


Case  2 


For  Case  2 

a)  d  =  0,  M  =  0 

b)  d  =  0.2, 0.5  cm 

•  Mf  =  0,0.1  W,0.2  W 

Mox  =  0 

•  M  F  =  .07  W  +  .07  Ba  +  .07  Au 
Mox  =  0 


c)  d  =  0.2,0.5cm 

Mf  =0  (100%  &  25%  density) 
Mqx  =  0.1  W,  0.2  W 


^/ence  Applications 
^  International  Corporation 

AnEiaphyM-OwimdCmfit/ty 


RR  91/797  h«A*a 
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CONCLUSIONS  OF  FSCATT  CALCULATIONS 


•  FSCATT  runs  show  10-30  percent  contrast  is  possible 
for  a  5  mm  jet  doped  with  high-Z  material,  or  for  a  5  mm 
jet  penetrating  a  doped  fluid 

•  Total  x-ray  transmitted  fluence  is  1  -  7  percent  of  incident 
fluence  depending  upon  doping  fraction,  geometry  and 
spectrum 

•  Film  response  data  indicated  marginal  exposure 
operation  at  150  kV,  acceptable  exposure  levels  at  300  kV 


Science  Applications 
Intamallonal  Corporation 

Ait  BuptofM-Owwrf  Comptny 


RR  91/797  haMa 


300  kV  EXPERIMENTS 


•  One  day  test  series  conducted  at  the  Alliant  Techsystems 
Proving  Ground,  Minneapolis 

-  HP  43733A  300  kV  FXR  system 

Kodak  XAR  film  using  TI-3  or  NDT-9  front  and  back 
intensifying  screens 

-  Automatic  and  hand  processing 

•  Two  mock-ups  used  as  x-ray  test  objects 

-  Water  filled  graphite  cylinder  with 
tungsten  carbide  doped,  gelled  water  on  axis 

-  Graphite  cylinder  filled  with  tungsten 
carbide  doped,  gelled  water  and  empty  glass 
tube  on  axis 


V 


Science  Af^lcatlon$ 
international  Corporation 


Rn91/B60h»/da 
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FLASH  X-RAY  RADIOGRAPHY  SUMMARY  AND  STATUS 


SUMMARY: 

•  FSCATT  x-ray  photon  transport  calculations  show  good 
contrast  in  tungsten  carbide  doped  SFI  or  CAP  geometries 

•  Carbon  composite  cyiinders  with  up  to  one  inch  wall 
thickness  may  be  used  as  the  x-ray  transparent  chamber 

•  300kV  FXR  experiments  confirm  code  predictions  and  serve 
as  initial  optimization  study 

•  PC  based  image  analysis  system  is  more  than  adequate  for 
the  task 

STATUS: 

•  Vendor  identified  for  composites;  FMC  working  on 
fabrication  specs  and  prototypes 

•  Burn  rate  chamber  to  be  ready  in  early  July,  1991 

•  Alliant  Techsystems  will  set  up  and  operate  the  FXR  system 


Scfence  A/^llcatlons 
tnt0mational  Corporation 


RR  91/860  htMa 
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X-RAY  TEST  OBJECT 


Tap  water 
20%  WC  doped 
gelled  water 


9  mm  OD,  7  mm  ID 
Glass  Rod 

•  Empty 

•  20%  WC  doped, 
gelled  water 


Graphite  Cylinder 
1.65"  ID 

•  3/4"  wall 

•  1/2"  wall 


Fill  and  vent  ports 


fT91/797tfa 

Science  Appllcetlona 
fntemetlonal  Corporation 

Afi  Ein)tey»*Or<iprf  Company 


RR  B1<797  hs>cla 


TEST  OBJECT: 

6RBPHITE  GV^LIKDER,  1.85"  ID. 
0.5"  I1RLLS.  GELLED  RATER  FILL, 
20^  HD  DOPED.  ID.  Smu  OD 

euptv  glass  rod  on  rkis. 

FILH:  KODAK  KAR 
FTSB:  60" 

OIFD:  3.75" 

F-SCREEN  TI-3 
B-SCREEH  TI-3 
300  kV  HP  FKR 
25  kV  CHARGE 
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TEST  OBJECT: 

ORRPHITE  CYLINDER,  1.85"  ID. 
0.5"  WRLLS.  GELLED  yfiTER  FILL. 
20^  HC  DOPED.  7mm  ID,  Smm  OD 
EHPTV  GLfiSS  ROD  ON  RKIS. 

Fflff:  KODAK  XAR 
FTSD;  80" 

OTFD;  3.75" 

F-SCREEK  KBT3 
B-SCREEH  ROTO 
300  kV  HP  FXR 
30  kV  CHARGE 


TEST  OBJECT: 

GRRPHITE’CVLINDER,  1.85^'  ID. 
0.5"  mtS,  GELLED  WATER  FfLL. 
20^  WC  DOPED.  7mffi  ID.  3mtn  OD 
EHPTV  GLASS  ROD  ON  RKfS. 

FILH;  KODAK  KARf 
FTSO;  48" 

OTFD:  3.75" 

F-SCREEU  11 -3 
B-SCREE8  TI-3 
300  kV  HP  FXR 
25  kV  CHARGE 
Z  HAND  PROCESS 
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TEST  OBJECT: 

SRfiPHiTE  CVLIHBER,  1.B5"  ID 
0.75"  HfiLLS.  iIRTER  FILL.  2oi  HC 
DOPED  OELLEO  WATER  IM  7mm  ID, 
Smm  OD  GLfiSS  ROD  ON  8X IS . 


FILH:  KODAK  XAR 
FTSD:  80" 

OTFD;  3.75" 
F-SCREEW  TI-3 
B-SCREEH  TI-3 
300  kV  HP  FXR 
25  kV  CHARGE 
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RETOUCHED  RADIOGRAPH  SHOWING 
SIMULATED  "PLASMA  JET" 


Science  Applications 
International  Corporation 

An  £mple)i**-Oitn9d  Cvmptny 


RRgi^g7hsAJa 
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f  EXAMPLE  OF  LOOK-UP  TABLE  OPERATOR 

GRADIENT  ENHANCEMENT 


RR  91/797  hsAfj 
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EXAMPLE  OF  HORIZONTAL 
GRADIENT  FILTER  APPLICATION 
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LASER  IN-BORE  POSITION  DIAGNOSTIC 


GOAL:  Provide  precise  projectiie  position  data  for  comparison 
to  interior  baiiistic  modei  predictions 

PROBLEM:  Microwave  radars  are  expensive  and  have  poor 

resoiution  at  low  velocity.  Optical  Doppler 
velocimeters  are  very  expensive  and  difficult  to 
field  in  ETC  gun  environments 

APPROACH:  Use  pulsed  laser  radar 

-  Inexpensive 

-  Robust 

-  Accurate;  reads  position  at  any  gun  velocity 

-  Easily  calibrated  m  situ 


Sefence  Appileationt 
Intarrmtlon^  CorporttUtm 

An  Employ— <fmt9dCoirpmf 


Rfl91/B60IW)te 
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PROJECTILE  VELOCITY  MEASURMENT  VIA 

LASER  RADAR 


LASER  RADAR  SYSTEM  PARAMETERS 


Laser  Wavelength 

820  nm 

Laser  Pulse  Width 

50  ps 

Collimated  Laser  Beam  Divergence 

0.15  mrad 

Absolute  Position  Accuracy 

2  cm 

Position  Jitter 

<  1 .5  cm 

Photodetector  Signal  Strength 

1.0  V 

Repetition  Rate 

40  kHz  - 1  MHz 
(A/D  board  options) 

$e/enc9  Applications 
International  Corporation 

Art  Cffipfcjw  Omtta  CwTpMy 


RR  91/860  hA/da 
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LASER  RADAR  COMPONENT  LIST 


PART 

MANUFACTURER 

COST 

150mW,  50ps,  1MHz  laser 
diode  @  820nm  with 

Hamamatsu 

$8500 

collimator 

35ps  photodetector 

Opto-Electronics 

3200 

Vibration  isolated 
optical  bench 

Newport 

2700 

Optical  components 
and  hardware 

Newport/Mel  les  Griot 

4900 

Analog  timing  electronics 

EG&G/ORTEC 

5250 

1Mhz,  16  Bit,  A/D  DAS  board 

Analogic 

3895 

Rack  mount  PC  system 

Cyber  Research 

3595 

DAS  software  (snapshot) 

Cyber  Research 

995 

Science  Applleetlone 
intematlenat  Corporation 

AnEir^hf-OmmlCan^tn/ 


$33,035  } 


RR  91/060  hft/d« 
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LASER  RADAR  SYSTEMS  ISSUES 


OPERATION  IN-BORE 

-  Projectile  blow-by  must  be  minimized 

-  Requires  precision  retro-reflector  mounted  in 
projectile  nose 

-  Requires  precision  opticai  alignment  for  each  shot 
TIMING  ELECTRONICS 

-  Millivolt  resolution  of  time-to-amplitude  converter 
signal  is  required 

-  Laser  system  and  timing  electronics  must  be 
electrically  isolated  and  heavily  shielded  against 
harsh  EMI  environment  of  the  ETC  gun 


V 


Science  App/tcalionB 
Intemetfonat  Coqyoration 

AnEnphyM 


HR  91/860  ha/d« 
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FIRST  PRINCIPLES  MODELING  OF  A  DMA  SOmm  ETC  GUN  DESIGN 


CC.  Hsiao,  G.  T.  Phillips,  and  F.Y.  Su 
Science  Applications  International  Corporation 
San  Diego,  CA  92121 


ABSTRACT 

The  SAIC  proprietary  model  BISON  is  a  first  principles  CFD  code  which  has 
recently  been  applied  to  the  study  of  the  interior  ballistics  phenomenology 
associated  with  a  class  of  ETC  gun  designs.  The  BISON  model  solves  the  NS 
equations  for  an  electrically  conductive  fluid  but  neglects  magnetic  forces. 
Bison  treats  multiple  phases,  turbulence,  multi-species,  and  chemistry  making 
it  well  suited  for  the  study  of  interior  ballistics. 

BISON  treats  the  advocation  terms  in  the  conservation  equations  using  a 
fourth  order  in  space,  second  order  in  time  FCT  algorithm.  Adaptive  gridding 
is  included  for  accuracy  and  economy  and  makes  use  of  both  continuous  and  discrete 
second  order  rezonlng  algorithms.  A  non-equilibrium  multi-phase  treatment  is 
available  in  BISON  through  a  discrete,  stochastic  particle  model  which  follows 
Lagrangian  computational  particles  which  are  fully  coupled  to  the  continuum  flow. 
In  addition,  the  multi-species  capability  in  BISON  allows  the  treatment  of 
equilibrium  two-phase  flow.  Multiple  fluid  species  are  treated  through  separate 
continuum  continuity  equations  for  each  species .  Real  equations  of  state  for 
each  species  are  provided  through  a  unique  non-iterative  multi-species  pressure 
and  temperature  solver.  Turbulence  is  treated  via  a  two-equation  (k-w)  model 
which  has  been  validated  for  a  variety  of  free  and  boundary  layer  shear  flows. 
Reactions  within  the  flow  are  computed  using  an  eddy-breakup  model  which  we 
validated  under  a  previous  DOE  IC  engine  program.  This  model  assumes  that  the 
chemical  kinetic  races  are  fast  compared  with  the  turbulent  mixing  rates  so  that, 
flame  propagation  is  mixing  controlled. 

This  presentation  will  demonstrate  the  effectiveness  of  the  BISON  model  in 
modeling  complex  interior  ballistics  processes.  A  test  case  will  be  shown  which 
describes  the  interior  ballistics  of  a  60  mm  ETC  gun  fixture  which  is  being  used 
in  the  DNA  ETC  Gun  Program.  Comparison  of  pressure  transducer  data  with  model 
predictions  clearly  shows  that  the  computational  approach  is  capable  of  capturing 
the  detailed  structure  present  in  the  experimental  data. 
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First  Principles  Modeling  of  a 
DMA  60mm  ETC  Gun  Design 


by 

CC  Hsiao 

Gary  Phillips  &  Fred  Su 

Science  Applications  International  Corporation 

July  9-11, 1991 
presented  at 

JANNAF  Workshop  ETC  Modeling  &  Diagnostics 
Aberdeen  Proving  Ground,  MD 


July  9.  1991 

•FMC  BISON  Model 


An  Smphyee 


High  Order  Numerics 


4th  Order  in  Space,  2nd  Order  in  Time 
FCT  Scheme 

High  Order  Adaptive  Gridding 
Arbitrary  Problem  Geometry  Specification 


\ 


First  Principle 
Physical  Models 

Plasma  Model 

Multi-Phase  Continuum  Model 
Multi-Phase  Discrete  Particle  Model 
Chemical  Reactive  Flow  Model 
Turbulence  Model 
Real  Equations-of-State 


Species  Transport  Equations 
Navier-Stokes  Equations 
Turbulence  Equations 


Governing  Equations 


Plasma  Model 
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July  9,  1991 


•FMC  Multi-Phase  Continuum  Model^igg 


An  Bfrfihy00  Own^d  C^mpsny 


Assumption:  Time  scale  to  reach  thermodynamic  equilibrium  is 
smaller  than  computational  time  step 


4th  order  FCT 


rnrnmr.mmmmm^m 


molecular  +  eddy 
diffusMties 


multi'Species  density  & 
energy  diffusion 


Chemical  reaction 


chemistry  & 
reaction  rate  model 
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^’^j^j.fXfjgi^Pgj^lf^lQ;^^  -1.- 


4o  Emfihy»9>^"fn*d  C<irrpMr\y 
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July  9,  1991 

•FMC  Turbulence  Model 


An  emptayt  Owfitd  CorrpAny 


Law  of  the  Wall 


Turbulent  Energy  &  Dissipation  Rate  Equations 


Advection  Production  Dissipation  Diffusion 

d(pK) .  d(pUiK)  (  ,  I  duk  d[,  , 

^1“  PS-P  w)ic-4  pK  -  (p+o  pe)^ 
dt  dxj  dx^  dxj  dXj 


apsJp+2o(-l1w|w2^ 

di  3xj  \dxj  /  J  I  dxj 


( p+OpE ) — - — 
dxj 


Eddy  Diffusivity 
Reynolds’  Stress  Tensor 
Turbulent  heat  flux  vectors 
Turbulent  species  flux  vectors 
Mixing  time  scale  for  reaction  rate  model 
Fluctuating  vel.  for  2-phase  particle  model 
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July  9,  1991 

f  Df  d  Real  Equations-of-State 


A  ri  gfnphy04  C^n4d  Ccmpdny 


Utilizing  existing  EOS  Packages 


DNA  EOS  Library 
Analytiq^  EOS 
Equilibrium  Chemical  Codes 
Bureau  of  Standard  EOS  data 


Verified  with  available 
experimental  data 


Range  of  Validity  covers  all 
.  possible  ETC  Gun  conditions 


High  density, high  pressure  liquid  statr 
Liquid-gas  phase  transition 
Supercritical  state 
High  energy  plasma  state 


Efficient  and  easy-to-use 
table  lookup  format 
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SPECIRC  VOLUME  (M^/KG) 


•  -  Data  from  Thermopyhsical  Properties  of  Matter, 
Purdue  University,  1982 
^  -  Critical  point,  from  Handbook  of  Chemistry 
and  Physics,  66th  Ed.  1986 


SPECIFIC  VOLUME 


o  -  Data  from  Hydrogen  Peroxide  Physical  Properties, 
FMC,  1969 

^  -  Critical  point.  Hydrogen  Peroxide  Physical  Properties, 
FMC,  1969 


®  -  data  calculated  from  NASA-CET89  cods 


W  ■  Critical  point,  from  Handbook  of  Chemistry 
and  Physics,  66th  Ed.  1986 


« -  Critical  point,  from  Handbook  of  Chernistry 
and  Physics,  66th  Ed.  1986 


2A6 


June  25.  1991 

•FMC 


Phenomenology  Experiments 

Burn  Rate  Tests 


An  Otvnatf  Company 


Plasma  Injector 
Closure  Plug 


•  Gain  an  understanding  of  plasma/fluid  interactions 

•  Allow  computation  of  gas  generation  rates 


•  Provide  data  input  for  model  verification 


June  25.  1  991 


Phenomenology  Experiments 

•FMC  Laser  Pyrolysis  Ignition  Tests 


An  Employe*  Owngg  Co/npnny 


Mirror ^ - 


Masl^ 


ZnSe  Window 


KCl  Window 


Oscilloscope 


Spin  Physics  2000 
High  Speed  Video 
Camera 


□ 


PC  Analysis 


Provide  fuel-oxidizer  chemical  kinetic  data 


Provide  ignition  delay  time  and  temperature  at  high  P 
Provide  data  for  chemical  reaction  computer  sub-model 
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June  25,  1991 


Phenomenology  Experiments 

Jet  Penetration  Tests 


An  £rtiQlOY««  Oim«d  Campan^r 


Pressure 

Regulator 


Mirror  ^ - 


ser  Sheet 


Spherical  Lens 
!|=>  Cylindrical  Lens 


Pressure  Transducers 


Gas  Suppf^: 
Reservoir 


Spin  Physics  2000  Pfojectiie 
High  Speed  Video 
Camera 


Measure  instantaneous  gas  jet  contours  and  entrainment  processes  at  near 
gun  pressures  (up  to  30  kpsi) 

Assess  shear  layer  stripping  and  droplet  generation 

Provide  input  data  for  model  verification 


249 


Jl  :  ■ 99 


fMC  60mm  Shot  26  jmr 


in  £rfphy00  Own^d  C^mpsny 


60  mm  Subscale  Gun  Fixture 


Electric  Energy: 

0.8  MJ 

Chemical  Energy: 

4.7  MJ 

Elec.  Pulse  Width: 

4  msec 

Projectile  Mass: 

2.75  kg 

Muzzle  Energy: 

1 .23  MJ 

moving  grid 


2D  axisymmetric  geometry 

SFI  annular  orifices  (3)  are  replaced  by  an 
annular  ring  having  the  same  surface  area 

5  species:  CgH^g,  H2O2,  H2O,  CO2,  and  Air 

are  included  in  the  computation 

4  material  regions:  plasma  channel, 
polyethylene  tube,  fuel  rich  (SFI),  and  fuel 
lean  propellant  (combustion  chamber),  are 
in  the  computation 
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July  9,  1991 

•FMC 


Capillary  &  SFI  Modeling 


An  Owntti  Coinptny 


Electrical  energy  data  from 
FMC  experiments  are  directly 
linked  to  the  computation 

Time  delay  for  fuse-wire 
breakdown  is  accounted  for 

Material  strength  of 
polyethylene  tube  is  ignored 

Plasma  channel  growth  is 
correlated  to  electrical  energy 
input 


1  2  3 

Time  (msec) 


HFliic  Chamber  Modeling 

_  y* ..  An  irrploymg  Owrwd  C^mptny 


Dbkjraior 


•  Flow  field  in  combustion  chamber  is  directly 
coupled  with  plasma  capillary  and  SFI 

•  Chemical  kinetic  rates  are  assumed  to  be 
fast  with  respect  to  turbulent  mixing  (implies 
combustion  is  mixing  controlled) 

•  Single  step  chemical  reaction  is  assumed 


HFiiic  Dynamics  Modeling 

An  tiTf^yt  Compu^i 


•  Dynamic  rezoning  follows  projectile 
motion 

•  Problem  geometry  Is  reconstructed 
according  to  new  grid 

•  Air  ahead  of  projectile  is  ignored  (  error 
<1%) 

•  Bore  friction  force  is  neglected 


Axial  Distance  (cm) 


.W  A.1. 


0.50  m3  LOD  ms  1.50  ms  1^.00  ms  2.50  ms  3,00  ms 
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Axial  Distance  (cm) 


Axial  Distance  (cm) 


0*50  ms  1-00  ms  1-50  ms  2«00  ms  2*50  ms  3,00  ms 
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0.50  ms  i.OO  ras  1.50  ms  S.OO  ms 
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Current  Design 


Design  Change  1 
Design  Change  2 


June  25.  1991 


'FMC 


SFI  Performance  Studies 


An  £mp/Dy««  Owned  Company 


Pressure  Transducer  P1 


0.5  1  1.5  2  2.5 

Time  (msec) 


[Change  1  Change  2 

50%  Orifice  Area  Reduction 


SFI  Face  Plate 


Current 

3 
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Current  Design 
Design  Ctiange  1 
Design  Change  2 


Current  Design 
Design  Change  1 


June  25. 1991 


SFI  Performance  Studies 


•FMC 


An  Employe  OHTWff  Cwnpjny 


Fuel  Jet  Density 


Density 


II » II 


Temperature 


Fuel  Jet  Temperature 


18000 
16000 
^  14000 
2  12000 
^  10000 


1.5  2 

Time  (msec) 


SFI  Face  Plate 


Current  iChange  1  Change  2 


50%  Orifice  Area  Reduction 


1  1.5  2 

Time  (msec) 
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Current  Design 


Design  Chartga  1 


Design  Chartge  2 


Current  Design 


Design  Change  1 


June  25.  1991 


>FMC  SFI  Performance  Studies 


An  Employe*  Onrrrstf  Company 


Jet  Velocity 


Mass  Injection  Rata 


Jet  Velocity 

Mass  Injection  Rate 


1  1.5  2 

Time  (msec) 

SFI  Face  Plate 


3  S  ^ 
a>  300 


Current  iChange  1  Change  2 


50%  Orifice  Area  Reduction 


1  1.5  : 

Time  (msec) 
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Current  Design 
Design  Change  1 
Design  Change  2 


Reaction  Cnamber 


Current  Design 
Design  Change  1 
Design  Change  2 


June  25,  1991 


*FMC 


SFI  Performance  Studies 


An  Qwrmd  Compuny 


Oxidizer  Injected 
into  Reaction  Chamber 


Fue 


Oxidizer  Reacts 


Oxidizer  Reacts 
in  Reaction  Chamber 


Fuel  Reacts 


in  SR 


in  Reecbon  Chamber 


Fuel  irijected 
into  Reaction  Chamber 


2  0.2 


Oxidizer 


nme  (msec) 


SFI  Face  Plate 


Current  •  Change  1  Change  2 


50%  Orifice  Area  Reduction 


Time  (msec) 


PHYSICS  OF  ETC  PIASHA-FLUID  INTERACTIONS 

B.  A.  Kashlwa,  H.  A.  Davis,  and  R.  J.  Tralnor,  Jr. 
Los  Alamos  National  Laboratory 
Los  Alamos,  NM  87545 


ABSTRACT 


Two  fundamental  elements  of  the  Electrothermal/Chemical  (ETC)  Gun  process 
are  discussed.  The  first  Involves  the  nature  of  discharge  arc  Itself,  and  the 
second  Involves  the  microphysical  Interaction  of  the  expanding  hot  gas  with  a 
cool  condensed  phase  propellant.  Magnetohydrodynamlc  calculations  suggest  the 
existence  of  a  so-called  filament  discharge  In  the  hot  expanding  gas.  These 
calculations  are  discussed.  Also,  recent  experimental  evidence  suggests  the 
Importance  of  Rayleigh -Taylor  Instability  in  the  development  of  a  cavity  in  a 
liquid,  at  prototypical  ETC  Gun  pressures.  The  latest  multidimensional  hydro- 
dynamic  analysis  of  this  Instability  effect  Is  presented. 
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B.  A.  Kashiwa,  H.  A.  Davis,  and  R.  J.  Trainor 
Los  Alamos  National  Laboratory 


PhyslGS  of  ETC  Plasma-Fluid  Interactions 


JANNAF  Workshop 
Jiily  9-11,  1991 
Aberdeen,  Md 
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Los  Alamos^  Goal  Is 
Predictive  Capability  For  ETC  Design 


(Log,  z^daoini®: 


Overview  of  Los  Alamos  ETC  Activities 

Working  Toward  A  Predictive  Capabiiity 


1 .  End-to-end  code  development 


2.  Benchmarking  expts.  on  20  mm  ET  gun 


3.  Plasma-working  fluid  interaction  expts. 


4.  Modeiing  of  P/W-F  experiments 


This  taik 


Hai  Davis 
next  taik 


5.  Projectiie  veiocity  diagnostics  — 


Dick  Bartsch 
tomorrow  AM 


Los  Alamos  ETC  Code  Uses 
Multidimensional  Hydrocode  Library 


•  End-to-end  code  combines  hydrocode  library  (CFDLIB), 
circuit  code  analysis  and  capillary  plasma  tube  model. 

►  Useful  engineering  analysis  tool 

►  Now  developing  into  wide-range  predictive  tool 

•  Gun  chamber  and  tube  treated  with  multidimensional  bode 
Which  treats  a  turbulent,  chemically-reactirig,  multisbecies 
fluid 

‘i* 

►  Predominantly  Eulerian  frame  used 

►  Two  propellant  burn  models  (Arrhenius  &  bulk  burn) 

,►  Can  access  variety  of  modules  (e.g.,  turbulence,  EOS) 

•  Reference  on  interfaces,  numerical  techniques,  models: 

B.  A.  Kashiwa  et  al.,  Proc.  27th  JANNAF  Combust.  Mtg.,  1990 


Los  Alamos  20  mm  ET  Gun 

Used  For  Code  Benchmarkina 


Voltacie 

4kV 

Energy 

■^"«^22-kJr''"'': 

CUrreriit 

10  kA 

lm|>eddhcie 

0.2  Q 

Power 

22  MW 

Max.  velocity 

0.6  km/s 

Rise  time 

0i6  ms 

Max.  pressure 

6  kpsi 

[L(q)i  /^Dii[ririi@. 


P-1/Dfi,HD 


Los  Alamos  ET  Gun  Diagnostics 


Pressure  Probe  No.  1  (kpsi) 

0.0  2.0  4.0  6.0 


Pressure  Probe  No.  2  (kpsi) 

0.0  3.0  6.0  9.0  12.0 


9zr 


Interferometer  Signal  (mV) 

-150.0  -100.0  -50.0  0.0 

I  I  .  ( 


OBSERVATIONS  AND  MODELING  OF 
FUNDAMENTAL  ELECTROTHERMAL  GUN  PHENOMENA 

H.  A.  Davis,  R.  R.  Bartsch,  B.  A.  Kaahlwa  and  N.  T.  Padlal 
Los  Alamos  Natioiial  Laboratory 
Los  Alamos,  NM  87545 


ABSTRACT 

We  have  begun  experimental  studies  coupled  with  computer  modeling  of 
fundamental  electrothermal  gun  processes.  In  Initial  work,  now  complete,  we 
used  a  20’mm-bore  steam  driven  electrothermal  gun  to  develop  diagnostics  and  to 
benchmark  Los  Alamos  ET  gun  computer  modeling.  Results  of  current  experiments, 
now  underway,  to  image  the  interaction  of  A  plasma  with  a  working  fluid  will  be 
reported.  The  object  is  to  determine  fluid  mixing  processes,  first  with  non¬ 
combustible  working  fluids  and  then  with  combi^stiole  fluids.  Initially,  the 
interaction  of  an  exploded  wire  plasma  with  water  is  being  studies  with 
time-resolved  X-ray  photography.  The  interaction  Is  contained  in  a  dielectric 
chamber  surrounded  by  a  steel  housing.  A  short  section  of  barrel  is  used  to 
allow  for  fluid  expansion.  First  shots  h^ve  been  fired  and  Images  have  been 
produced. 
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H.A.  Davis,  R.R.  Bartsch, 
B.A.  Kashiwa,  and  N.T.  Radial 
Los  Alamos  National  Laboratory 


•  ■> 


of  Fundamental  Electrothermal 


JANNAF  Workshop 
9-11  July,  1991 
Aberdeen,  Md 


Los  Alamos  Research  Goal 

Investigate  Fundamental  ETC  Phenomena 


•  Understand  Mixing  and  Control 

Through  Well  Diagnosed 
Experiments  (Imaged  Chambers) 

s  •  Develop  Predictive  Computer 

Models  Using  Experimental  Observations 

•  Approach  -  Begin  with  Plasma  Embedded 

in  Working  Fluid  (Inert  and  Combustible) 


Electrodes 


Projectile 


Insuietor 


Liquid  Filled  _ 
Transparent  Chamber 


Wirrdow 


Steel  Housing 


P-1/0B.HD-2« 
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P-1/OB-17 


ELECTROTHERMAL  GUN 
WITH  TRANSPARENT  CHAMBER 


Stainless  Steel 


Fiber 


Lexan 

Chamber 

(liquid  filled) 


\  Lexan 
y  Bullet 


Spoked 

Electrode 


1 
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Electrical  Drive  Circuit 


+ 


2.74 

mF 


Si 


100  ^iH  18  mQ 


ETGun 

Load 


Load  Voltage 
Load  Current 
Energy  Delivered 
Peak  Pressure 
Projectile  Velocity 


4-8  kV 
6  kA 

10  kJ  in  0.5  ms 
10-20  kpsi 
250-300  m/s 


■\  f 


P-1/DB.HD-27 


Plasma/Working-Fluid  Interacterion  Experiment 


LORD  UOLTRGE  CKUOLTS> 


Load  Voltage  Waveform 


CURREHT  CKniiPS> 
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CHRMBER  PRESSURE  CKPSI> 


Radiogi  aphic  Configuration 

Copper  Fiber-0.33  mm  Diameter 


Orientation 


Time 


2.5  cm 


Projectile 


4.0  cm  — 


SHOT  20 
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What  is  the  Mechanism  Driving 
A  Large  Scale  Instability? 

•  Simulate  discharge  using  simple  heat  source 

/  . 

to  replace  plasma  model: 


ep  =  A  e“''  [l  +  B  sin^ 


radial  axial 

deposition  deposition 


Time-dependant 

modulation 


•  Parameters  chosen  to  match  experimervt.  Amplitude 

chosen  to  give  total  power  deposition  of  20  k J/ms . 

[L®§  /^DiiDTfa® 


RADIUS 


Perturbation  Is  Stabilized 


t  =  0.1  ms 


t  =  0.2  ms 


Z-AXIS 


Z-AX)S 


Chamber  Pressure  (kpsi) 


Chamber  Pressure  vs.  Radius 
at  Different  Times 


Radius  (i  =  10)  [cm] 
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Chamber  Pressure  vs.  Time 


Perturbation  Is  Unstable  With  Time  Modulation 


Perturbation  Is  Unstable  With  Time  Modulation 


0.18  ms 


0.20  ms 


Chamber  Pressure  (kpsi) 


Chamber  Pressure  vs.  Radius 
At  Different  Times 


Without  Modulation  With  Modulation 


0.0  0.2  0.4  0.6  0.8  1.0  1  2  0.0  0.2  0.4  0.6  0.8  1  0  1.2 


Radius  (i  =  10)  Icml 


Radius  (i  =  10)  Icml 


f BARREL  PRESSURE 


Chamber  Pressure  vs.  Time 


Summary 


Pressure  and  Electrical  Measurements 
Suggest  Pressure  Waves  are  Important 

Radiographs  Show  Instability 

(axial  wavelength  about  equal  to  radius) 

Calculations  Show  Damping  of  Purely 
Spatial  Perturbation 

Addition  of  Power  Modulation  in  Calculations 
(at  Chamber  Transit  Period)  Produces  Instability 
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Future  Directions  For  ETC 
Computational  Developmnet 


•  Incorporate  detailed  plasma  dynamics 

►  Sensitivity  of  conductivity  to  P  and  T 

►  Can  P  waves  be  sustained  by  sensitivity  of  ohmic 
heating  to  pressure  fluctuations? 

•  Couple  circuit  dynamics 

•  Examine  potential  of  driving  instabilities  in  propellants 

•  How  can  we  avoid  resonance  between  pressure  waves 
and  plasma  dynamics?  Begin  examining  different  kinds 
of  electrode  configurations 


Future  Experiments 


•  Study  Different  Configurations 
y  Vary  Geometry 

y  Combustible  Working  Fluids 
y  Vary  Power  Input 

•  Improve  Measurements 

y  Orthogonal  X-Ray  Views 
y  Time  Sequence  on  Single  Shot 
y  Time  Resolved  Photography 
y  Time  and  Spatially  Resolved 
Spectroscopy 


Intentionally  left  blank. 
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ELECTROTHERMAL- CHEHICAL  GUN  MODELING 


David  W.  King,  Olin  Rocket  Research 
and 

Hugh  McElroy,  Olin  Ordnance 


ABSTRACT 

A  presentation  will  be  made  describing  the  nature  of  ET-C  Gun  Modeling. 
The  key  physics  regarding  combustion  and  arc  physics  will  be  discussed  with 
regard  to  implementing  a  2-D  CFD  model  of  the  gun.  At  OLIN,  both  liquid  and 
solid  propellant  2-D  modelling  technology  is  being  developed.  The  status  of 
liquid  propellant  model  will  be  presented  along  with  example  calculations. 
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Olin  ROCKET  RESEARCH  COMPANY 

DtfiHse  SYsreus  Cfioup 


ELECTROTHERMAL-CHEMICAL  GUN 
MODELING 


ELECTROTHERMAL-CHEMICAL  GUN 

MODELING 

HUGH  A.  McEIROY 
OLIN  ORDINANCE 

DAVID  Q.  KING 

ROCKET  RESEARCH  COMPANY 
OLIN  DEFENSE  SYSTEMS  GROUP 

JULY  11, 1991 
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Olin  ROCKET  RESEARCH  COMPANY 

Dtrenu  srsrcMi  caoua 


ELECTROTHERMAL-CHEMICAL  GUN 
MODELING 


ET-C  CODES  AND  APPLICATIONS 


EQUATION  OF  STATE  FOR  HIGH 
PRESSURE  CHEMICAL 
EQUILIBRIUM 


HUNTER  &  BLAKE 
give  Impetus,  y 


BKW  gives  Curve  Fit  yiP,T) 


e  1  . 


t-ii 

Ul 


-S:';ar  s!*"  I 


^2-0  =- 
---  -L  ■“  “  s  Kiss 


USED  FOR  INTIAL  PROPELLANT 
SCREENING 

PROVIDES  IMPETUS  AND  FLAME 
TEMPERATURE 


*  MANY  A r’pr::;'x.M ATfONa.  ■ : 

-  - ^  -..v"':  h: : 

-  ilI 'O ^ - o-«  ^ 


*  W.VA  '1  f;t 


f  L.  I.I-N  ^  r  E"  L^lr-:=CE^f1Hi-’Ai=  i 


*  sTh  '  =  :  IL;7j[-:\'-;^'WT-rI 

*  y  sT  1  ■■■  =-  i'=  J I  =  ■  J-= !  L 
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Olin  ROCKET  RESEARCH  COMPANY  ELECTROTHERMAL-CHEMICAL  GUN 

DtF£H%t  SrSTEyS  GAOI/P  MODELING 


KEY  PHYSICS  ISSUES  IN  ETC  MODELING 


•  1D,  QUASH  D,  2D 

-  1 D  MODELS  OR  QUAS1 1 D  REGIONS  HAVE  INFINITE  SOUND  SPEED  IN  RADIAL 
DIRECTION,  FINITE  IN  AXIAL  DIRECTION 

-  PECULIAR  EFFECTS  ARE  OBTAINED  IN  RADIAL  DIRECTION  AT  INTERFACE  OF 
ID  TO  2D  REGION 

•  2D  MODEL  CAN  REPRESENT  RADIAL  VARIATION,  WHICH  PLAY  A 
ROLE  IN  HIGH  FREQUENCY  PRESSURE  WAVES 

U) 

Ch 
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Clin  ROCKET  RESEARCH  COMPANY  ELECTROTHERMAL-CHEMICAL  GUN 

OertNSe  SrsTtMS  G^OUP  MODELING 

EQUATION  OF  STATE _ _ 

•  REAL  GAS  IN  EQUILIBRIUM 

•  MULTIPHASE 


DQK;3 
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Olio  ROCKET  RESEARCH  COMPANY  ELECTROTHERMAL-CHEMICAL  GUN 

otftNse  srsrtkis  CKOU/f  MODELING 


COMBUSTION _ _ 

•  GAS  PHASE  IN  CHEMICAL  EQUILIBRIUM 

-  KINETIC  RATES  VERY  HIGH  AT  HIGH  PRESSURE 

•  EVOLUTION  OF  SOLID  OR  LIQUID  TO  GAS  PHASE  CONTROLS 
COMBUSTION  RATE 

-  DETAILED  PROCESS  INVOLVES  DIFFUSION  OF  HEAT  AND  SPECIES 

-  CAN  BE  CHARACTERIZED  BY  A  TIME  OR  LENGTH  SCALE  TYPICALLY  MUCH 
SMALLER  THAN  DIMENSIONS  OF  COMPUTATIONAL  CONTROL  VOLUME 

-  SURFACE  AREA  DEPENDENT 

•  REPRESENTATIONS  IN  A  COMPUTATIONAL  CONTROL  VOLUME 

-  VOLUME  AT  P,T  CONTAINS  UNBURNED  AND  BURNED  PROPELLANT 

-  BURN  RATE  GIVEN  BY  EXPERIMENTAL  PRESSURE  OR  TEMPERATURE 
CRITERION 

-  MASS  AND  HEAT  TRANSPORT  ACROSS  CONTROL  VOLUME  SURFACES 
STRONGLY  AFFECT  COMBUSTION  INSIDE  CONTROL  VOLUME 

•  MASS  AND  ENERGY  TRANSPORT  GOVERNED  BY  MIXING 
LENGTH  TURBULENCE  MODEL 

DQK:4 


Olin  ROCKET  RESEARCH  COMPANY  ELECTROTHERMAL-CHEMICAL  GUN 

biFtttse  SYsreus  cboup  MODELING 

ARC  HEAT  INPUT _ _ _ _ 

.  KILOAMP  DISCHARGE  AT  KILOBAR  PRESSURES  HAVE  VERY 
SMALL  RADIUS 

-  RADIUS  GIVEN  BY  PRESSURE  AND  ENERGY  BALANCE  ACROSS  ARC  TO 
FLUID  BOUNDARY 

-  TYPICAL  R  <0.5  mm 

•  VERY  SMALL  VOLUME  OF  FLUID  IS  IONIZED 

•  INTENSE  RADIATION  FROM  ARC  IS  TRAPPED  IN  VERY  SMALL 
DISTANCE 

.  IONIZED  FLUID  RECOMBINES  RAPIDLY  AT  KILOBAR  PRESSURES 

•  EXPECT  VERY  LITTLE  EFFECT  ON  CHEMISTRY  FROM  IONS 

•  TREAT  ARC  AS  LOCALIZED  HEAT  SOURCE 
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Olin  ROCKET  RESEARCH  COMPANY  ELECTROTHERMAL-CHEMICAL  GUN 

oiUMse  trsrtAis  gmup  MODELING 

STATUS _ _ _ 

•  APPROACH 

-  2D  CODES  DEVELOPED  BY  LOS  ALAMOS  NATIONAL  LABORATORIES  UNDER 
CONTRACT  TO  ROCKET  RESEARCH  COMPANY/OLIN  DEFENSE  SYSTEMS 
GROUP 

-  IN  HOUSE  MODIFICATION  FOR  SPECIFIC  APPLICATIONS 

.  LIQUID  PROPELLANT  2D  MODEL 

-  OPERATIONAL 

-  MULTIPHASE,  SINGLE  VELOCITY  &  TEMPERATURE 

-  VAUIDATION  WITH  EXOTHERMIC  PROPELLANTS  UNDERWAY 
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Oliii  ROCKET  RESEARCH  COMPANY  ELECTROTHERMAL-CHEMICAL  GUN 

otfettse  srstuts  gkoup  MODELING 


STATUS  (CONTINUED)  _ _ .-  .y;.:: 

•  SOLID  PROPELLANT  2D  MODEL 

-  READY  FOR  INITIAL  TESTING 

-  MULTIFLUID  -  SEPARATE  TEMPERATURES  &  VELOCITIES  FOR  GAS  AND 
SOLID  FLUIDS 

-  HANDLES  SOLID  PROPELLANT  AND  FRACTURED  PIECES  OR  "RUBBLE"  WITH 
CLOSE  RACKED  BED 


N» 
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Olin  ROCKET  RESEARCH  COMPANY 
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ELECTROTHERMAL-CHEMICAL  GUN 
MODELING 

EXAMPLE  CALCULATION  -  STEAM  GUN  WITH  CENTERLINE  HEATING 

DENSITY  CONTOURS 


t4).58  n» 


T-6.3  « 


T-0.2  m 


T»0.3  ms 

<l - 
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Olin  ROCKET  RESEARCH  COMPANY  ELECTROTHERMAL-CHEMICAL  GUN 

D£r£use  SYsnus  ottoup  MODELING 

EXAMPLE  CALCULATION  -  STEAM  GUN  WITH  CENTERLINE  HEATING 

PRESSURE  CONTOURS 


PRESSURE  TAPS 


(cm) 
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DtFiNse  iYsreus  gaoup  MODELING 


APPLICATIONS _ 

•  GAS  GENERATORS 

•  ET-C  MODELING 

-  VARtOUS  GEOMETRIES 
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ELECTROTHERMAL-CHEMICAL  GUH 
MODELING 


CONCLUSIONS 


•  LIQUID  &  SOLID  PROPELLANT  MODELS  ARE  OPERATIONAL  AND 
VALIDATION  IS  UNDERWAY 

.  APPROACH  YEILDS  A  CAPABILITY  THAT  IS  EASILY  ADAPTED  TO 
PROBLEM  SPECIFIC  SPECIALIZATIONS 
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RAILGUN  RESEARCH  RELEVANT  TO  ELECTROTHERMAL  GUNS 


Jad  H.  Batteh 

Science  Applications  International  Corporation 
Marietta,  GA  30062 


ABSTRACT 

Some  aspects  of  the  modeling  and  diagnostics  of  electrothermal  guns  are 
common  to  railguns,  as  well.  The  objective  of  this  presentation  is  to  compare 
and  contrast  the  two  technologies  and  to  provide  a  reference  guide  to  the  recent 
relevant  literature  describing  rallgun  research.  Particular  emphasis  is  placed 
on  the  following  areas  -  the  calculation  of  thermodynamic  and  transport  properties 
in  very  dense,  relatively  low  temperature  plasma;  on-board  projectile  diagnostics; 
spectroscopic  analysis;  and  the  confinement  and  stability  of  high-pressure  arc 
discharges. 
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Railgun  Research  Relevant  to 
Electrothermal  Guns 


Jad  H.  Batteh 

Science  Applications  International  Corporation 
1519  Johnson  Ferry  Road 
Suite  300 

Marietta,  Georgia  30062 
(404)973-8935 

Presented  at  JANNAF  Workshop  on 
ETC  Modeling  and  Diagnostics 
July  9 -11, 1991 
Aberdeen  Proving  Ground,  MD 


Objective 


•  Review  railgun  research  relevant  to  ETC 
technology 

-  Analysis  of  plasma  jets 

-  Diagnostics 

•  Provide  references  to  the  relevant  literature 


^4/L. 
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Basic  References  for  Railgun  Physics 


•  Proceedings  of  the  Symposia  on  Electromagnetic  Launcher 
Technology  published  by  IEEE  Transactions  on  Magnetics 

1.  Vol.  18,  January  1982 

2. Voi.  20,  March  1984 

3.  Vol.  22,  November  1986 

4.  Vol.  25,  January  1989 

5.  Vol.  27,  January  1991 

•  Special  Issue  on  Electromagnetic  Launchers,  iEEE  Transactions  on 
Plasma  Science.  Vol.  17,  June  1989. 

•  “Proceedings  of  the  Electromagnetic  Launcher  Workshop  II  - 
Diagnostic  Techniques”,  AFATL-TR-87-65,  January  1988. 


Comparison  of  Plasma  Flows 


*o_  thermodynamic  pressure 
^  ~  magnetic  pressure 
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MOPEUWG  DIFFICULTiES 


INCOMPLETE  UNDERSTANDING 
LIMITS  OUR  ABILITY  TO  CONTROL 
ARMATURE  PERFORMANCE 


•  MECHANISMS  WHICH  ARE  CONTROL  ARC  MASS 
AND  COMPOSITION  ARE  NOT  WELL  UNDERSTOOD 

-  COMPLEX  INITIATION  PHYSICS 

-  ARMATURE  EXCHANGES  MATERIAL  WITH 
BORE  THROUGH  ABLATION  AND  PLATING 

•  CONVENTIONAL  ELECTRIC  AND  THERMAL 
TRANSPORT  MODELS  ARE  INADEQUATE 
(NONIOEAU.  PLASMA) 

•  PLASMA  REPRESENTS  A  NON  -  STEADY,  3-D. 
TURBULENT,  POSSIBLY  MULTI  -  PHASE  MHO  FLOW 

•  PLASMA  HAY  BE  SUSCEPTIBLE  TO  INSTABIUTIES 


DIAGNOSTIC  DIFFICULTIES 

•  LIMITED  DIAGNOSTICS  CURRENTLY  AVAILABLE 

•  RELIABLE,  HIGH  FlDELfTY,  CONTINUOUS 
MEASUREMENTS  OF  ACCELERATION  ARE 
UNAVAILABLE 

•  HIGH  PRESSURES  IMPEDE  OPTICAL  DIAGNOSTICS 
OF  PLASMA  CORE 


•  PRESENCE  OF  STRONG  TRANSIENT  HELDS 


SA90160J>6 


PLASMA  PARAMETERS 


NONibEAL  PARAMCTER 


_  Colomb  Energy 
Thermal  Energy 


£>■(1+2)  n 


4n£  kl  l  (i+a)2/3  j 


NUMBER  OF  PARTICLES  IN  A  DEBYE  SPHERE 


Np  «  (f) 


DEGENERACY  TEMPERATURE  RATIO 


(^) 


rtie  k 


TYPES  OF  PLASMAS 


REQUIREMENT 

IMPLICATION 

IDEAL 

Y  «1 

Neglect  Coulomb  energy 
relative  to  thermal 
energy  in  plasma 

NONIDEAL 

Y  ^  0.2 

Must  account  for  Coulomb 
energy  in  determining 
properties 

QUANTUM 

Free  electrons  no  longer 
described  by  Boltzmann 
statistics,  and  plasma 
properties  must  be  based  on 
Fermi— Dirac  quantum 
statistics 

!»ry 


Scitnc*  ApplletUon*  lnl»fnMllon»l  Corpor»ll«n 


Electron  Number  Density 


METAL-PLASMA  ARMATURES  ARE 
TYPICALLY  NONIDEAL 


Temperature  -v  K 


Assumes  Singly  Ionized  Plasma 


Se<«>ie»  AfipHcatlon*  Cprpoftllon 
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IMPACT  OF  NONIDEAL  PLASMA  BEHAVIOR 


•  SIGNIFICANT  DECREASE  IN  THE  PLASMA  CONDUCTIVITY  BELOW 
THAT  PREDICTED  BY  SPITZER'S  FORMULA 

•  LOWERING  OF  THE  IONIZATION  POTENTIAL 

•  DECREASE  IN  GAS  PRESSURE  WHEN  COMPARED  TO  THE  IDEAL 
GAS  EQUATION  OF  STATE 


Selected  References  for 
Nonideal  Plasma  Properties 


*  “Electrical  Conductivity  of  Nonideal  Quasi  -  Metallic  Plasmas,”  A.H. 

Khalfaoui,  IEEE  Trans,  on  Plasma  Science  PS-12, 179  (1984) 

•  “The  Electrical  Conductivity  of  Plasma  in  a  Wide  Range  of  Charge  Densities,” 
Yu.  K.  Kuriienkov  and  A. A.  Vaiuev,  Beitr.  Plasma  phys.  24, 161  (1984). 

•  “Physical  Properties  of  High-Pressure  Plasmas  of  Hydrogen  and  Copper  in 
the  Temperature  Regime  5000  -  60000  K,  “P.  Kovitya,  IEEE  Trans,  on  Plasma 
Science  PS-r5, 587  (1985). 

*  “Electricai  Conductivity  of  Nonideal  Plasmas,”  R.J.  Zollweg  and  R.W. 
Liebermann,  J.  Appl.  Phys.  62,  3621  (1987). 


“Thermpdynamic  and  Electrical  Properties  of  Railgun  Plasma  Armatures,” 
(^^#lii[|^®r  and  J.K  E^eh,  IEEE  Trans,  on  Plasma  Science  17, 439  (1989) 


it  C|MlttilPn|iiand  Thermodynamic  Functions  of  Weekly  Nonideal 
Plasma,”  R.B.  Moh^and  J.G.  Gilligan,  J.  Appl.  Phys.  68,  5044  (1990). 


1/R  (theoretical) 


Comparison  of  Measured  and  Predicted 
Conductance  for  Plasma  Capillaries 

•  Based  on  quasi-steady  model  of  Powell  (in  Technology  Efforts 
in  ET  Gun  Propulsion,  Vol.  2, 1989) 

•  Predictions  assume  choked  flow 

•  Measured  values  denote  conductance  in  units  of  mhos 

•  Non-ideal  model  based  on  Kurilenkov  and  Valuev 

Uj  •  .  r  . 

06 


Ideal  Non-Ideal 
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Factors  That  Influence  Current  Distributions 
in  Confined  and  Unconfined  Discharges 


•  Solid  boundaries 

-  raiigun  bore 

-  capillary  wall 

•  Dynamic  (pressure)  equilibrium 

-  discharge  into  liquid 

•  Electric  field  distribution 

•  Conductivity  distribution 

•  Flowfield 

J  =  o  (E  +  V  X  B) 

•  Magnetic  confinement 

-p<1 


SflJL. 


Instabilities 


MUZZLE  VOLTAGE 
tARMATURE^FRClENCY) 


□ 


I  ~  I  FLASH  X-RAYS  & 


MAGNETIC  "G-DOT"  COILS 
(SENSE  ARMATURE  PASSAGE 
&  CURRENT  DISTRIBUnON 
IN  ARMATURE)  ^  ^ 


BREAK  SCREENS 
(PROJECTILE  VELCiClTY 
&  INTEGRITY) 


MEINEL 

SPECTROGRAPH 

(COMPO^ION,  TEMP 
&  DENSITY  OF  MUZZUE) 


I  PRESSURE  GAGES 
(IN-BORE  PRESSURES) 


BREECH  CURRENT 
(ROGOWSKI  COIL)" 


BREECH  VOLTAGE 
(LAUNCHER  EFFICIENCY) 


CONDUCTIVITY  PROBES 
(ARC  LENGTH) 


SPECTROGRAPH/ 

POLYCHROMATOR 

(COMPOSmON.  TEMP 
&  DENSnY  OF  ARMATURE) 


THE  GOAL  OF  THE  DIAGNOSTIC  DEVELOPMENT  EFFORT  IS  TO  ACQUIRE 
THE  NECESSARY  RAW  INFORMATION  THAT  CAN  BE  USED  TO  DETERMINE 
COMPONENT  CHARACTERISnOS  AND  EML  PERFORMANCE 


SA80161B-18 
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Overview 


Diagnostics 

Spectroscopy 
of  plasma 
armatures 


X-ray 

absorption  of 
plasma  armature 


of  Advanced  Railgun  Diagnostics 


Measured  or 
inferred  Parameters 

•  Composition 

•  Temperature 

•  Flow  Structure 


•  Density 

•  Composition 


Issues 

•  Short  optical  depths 

•  Data  interpretation 
transient 

•  Coating/shielding 


•  Initial  estimate  of 
components 


Selected 

References 

Keefer,  et  al. 
-Ref.  1,p217 

-  Ref.  2,  p.  295 
•  Ref.  3,  p.  360 

Clothiaux,  et  al. 
-Ref.  1,p.  199 

Clothiaux,  et  al. 

-  Ref.  1,  p.  139 
-Ref.  1,p.  165 
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Overview  of  Advanced  Railgun  Diagnostics 


'Measured  or 

-  '  •  ;  -I.  ; 

rSelected 

Diagnostics 

inferred  Parameters 

Issues 

References 

Magnetic 

•  Armature  current  - 

•  3-D  structure 

Cobb,  et  al. 

B-dots  ' 

distribution 

•  Qu'asi-steady 

-  Ref.  1,  p.  189 

•  Armature  iocation 

assumptions 

-  Ref.  3,  p.  507 

•  Orientation  ‘ 

1':.; ;  ^  ^ 

"Jamison,  et  al. 

-  Ref.  2,  p.  256 
Parker 

»  • 

r’s:. 

-  Ref.  3,  p.  487 
Smith,  et  al. 

-Ref.  3,  p.  501 
Bouvier 

-  Ref  .  3,  p.  516 

On-board 

•  Projectiie 

•  Soft-catch 

Fernandez,  et  al 

diagnostics 

acceleration  and 

•  Vibrations 

-Ref.  1,  p.  185 

velocity  profiles 

•  Mass/volume  of 

Littrell,  et  al. 

package 

■  Ref  4 

Grzesik,  et  al. 

-  Ref.  1,  p.  147 

Canpmnf 
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Overview  of  Advanced  Railgun  Diagnostics 


Measured  or 

Diagnostics  inferred  Parameters 
Visar  *  Veiocity  profiies 


Caiorimetry  •  Total  thermai 

energy  deiivered 
to  surtace 


Issues 


Selected 

References 


•  Shock  heating  of  Asay,  etal. 

ambient  gas  •  R®*-  2,  P-  46 

•  Vaporization  of 
projectiie/bore 

•  Blowby 

•  Projectile  deflectiori/ 
deformation 

•  Mechanical  vibrations 


•  High  powers 


Derbidge,  et  al. 
•Ref.1,p.202 


Radiometry  ‘Heat  flux 


•  High  powers 

•  Vibration  sensitivity 


Ibrahim 
-  Ref  5,  p.  2045 


Reference  List  for  Diagnostics  Table 


1.  IEEE  Trans,  on  Magnetics  g7,  Jan.  1991 

2.  IEEE  Trans,  on  Magnetics  25,  Jan.  1989 

3.  IEEE  Trans,  on  Plasma  Science,  June  1989 

4.  Third  European  Symposium  on  EML  Technology 

5.  J.  Phys.  D:  Appl.  Phys.  ^  (1980) 
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Summary 


•  Some  railgun  research  is  relevant  to 
ETC  technology 

-  Plasma  physics 

-  Plasma/soiid/liquid  interactions 

-  Diagnostics 

-  Data  analysis 

•  Railgun  literature  is  well-documented  and 
readily  available 

»  Expect  surprises 


^nir 

OVMtri 


INTENTIONALLY  LEFT  BLANK. 
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IN-BORE  POSITION  AND  VEtX)CITY 
MEASUREMENT  TECHNIQUES 

R.  Richard  Bartsch  and  Harold  A.  Davis 
Los  Alamos  National  Laboratory 
Los  Alamos,  NM  87544 


ABSTRACT 

A  mode-launcher  has  been  developed  which  permits  in-bore  microwave  inter¬ 
ferometry  in  the  lowest  propagating  mode  (TE^i) .  This  technique  eliminates 
multipath  and  multiple  mode  problems  which  can  occur  with  a  launcher  external 
to  the  gun  barrel,  and  may  result  in  position  data  which  is  clean  enough  to  look 
at  details  of  the  acceleration  of  the  projectile.  Tests  with  a  20inm  steam  gun 
will  be  discussed. 


This  work  performed  under  the  auspices  of  the  U.S.  Department  of  Energy. 
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and  Harold  A.  Davis 


Los  Alamos  National  Laboratory 

In-Bore  Position  and  Velocity 
Measurement  Techniques 


JANNAF  Workshop 
9-11  July,  1991 
Aberdeen,  Md 


Outline 


•  Fundamental  Mode  Microwave 
Interferometer  (20  mm  Test  Gun) 
y  Developed  and  Demonstrated 
y  Position  Measurement 

•  R.F.  Modulated  Laser  Radar 

y  Developed  and  Demonstrated 
@  500  MHz 

y  Designed  1  GHz  System 
y  Position  Measurement 

•  VISAR  (Velocity  Interferometer 
System  for  Any  Reflector) 

y  Developed  at  SNL 
y  Refined  and  Extended  at  LANL 
y  Velocity  Measurement 


CONVENTIONAL  EXTERNAL  LAUNCHER 


a)  INEFFICIEKT  d)  EASY  TO  SETUP 

b)  MULTIPLE  PATHS 

c)  MULTIPLE  MODES 


PROJECTILE 


GUN  BARREL 


RECTANGULAR 

WAVEGUIDE 


INPUT/OUTPUT 


FUNDAMENTAL  MODE  LAUNCHER 

a)  MECHANICAL  DESIGN  b)  SINGLE  MODE 
TO  WITHSTAND  BLAST  c)  EFFICIENT 
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P-1/DB-25 


f  =  Frequency 

k  s  Propagation  constant  =  2;c/guide  wavelength 
fop  =  Operating  frequency  of  conventional  interferometer 
fop  s  Operating  frequency  of  fundamental-mode 
interferometer 

PROPAGATION  OF  WAVEGUIDES  MODES 

L@§ 
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92-9a/l^-d 
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0’09-  O'OOl* 

(AUl)  leuOjS  J9l9UJOJ8^J8)U| 


lONAL  <niLLlU0LT8> 


Microwave  Interferometer  Signal 


1 


SINEWAVE  MODULATION 
(CW  source) 


PROJECTILE 


GUN  BARREL 


REPLACEABLE 

MIRROR 


MODULATED 
LASER  BEAM 


LASER  RADAR 
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PROJECTILE 


fsa 


DIRECTIONAL 

COUPLER 


rf  OSC. 
-  1  GHz 


OUTPUT 

SIGNALS 


AMP 


STEPDOWN 

XFORMER 


m 


4  AB  cos  0 


4  AB  sin  9 


l-l  PHASE 
DISCRIMINATOR 
(correlator) 


B 


1" 

LASER 

-850  nm 

AM  LASER 

L 

DRIVER 

LASER 

OUTPUT 

FIBER 

OPTIC 


DIODE 

BIAS 


REPLACEABLE 
REFLECTOR 


DETECTOR 
DIODE 


OUTPUT 
LENS 


Projectile 

pIsj^ieioBinerit 


|09V»loi^  6t  500IHHz%tl»m| 
0^  jL^h^>:p4:UNL 


»f  MODUI.ATFI>-P.ASER-BEAM  POSITION  SENSOR 


P-1/DB-C23 


PROJECTILE 


GUN  BARREL 


REPLACEABLE 

MIRROR 


IMAGING  DELAY 
SYSTEM  (>-1.6  ns) 


DIFFUSE  REFLECTOR 
(non-optlcal  quality 
surfacs) 


OVERLAY  OF  IMAGE 
AND  DELAYED  IMAGE 


CW 

INPUT 

BEAM 


PM  TUBE 


V9SAR 


FRINGE  COUNT 
GIVES  VELOCITY 
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Lhe  VISAR  Vgloclijj  Hi434 


DIFFUSELY  REFLECTING  SURFACE 


PROJECTILE 


GUN  BARREL 


CnaL 


PHOTOMULTIPLIERS 


POLARIZER 


BEAM 
A  SPUTTER 


POLARIZING 
BEAM  SPLITTER 


Ml  large  beam 
SPLITTER 


1/4  WAVE 
PLATE 


ETALON 


SCHE-ViAT’C  :F  a  -'iSAR  i 


ref.  Barker  and  Hollenbach, 

J.  Appl.  Phys,  Vol.  43,  No.11,  4669 
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220-qa/t-d 


0  4  i::-8  12  16  20  24 

TIME  (ms) 

FRINGE  DATA  FOR  PROJECTILE  MOTION 

ref.  Barker  and  Hollenbach, 

J.  Appl.  Phya.  Vol.  43,  No.1 1 , 4669  L©®  /?\0il(nnl© 
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Summary 


•  Fundamental-Mode  Microwave 
Interferometer 

•  R.F.  Modulated  Laser  Radar 

•  VISAR 


INT^TIONALLY  left  BLAiW. 
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IN- BORE  ACCELERATION  HEASDREMENTS  WITH  AN 
INSTRUMENTED  RAILGUN  PROJECTILE 

Donald  M.  Llttrell,  Keith  A.  Jamison  and  Glenn  E.  Rolader 
U.S.  Air  Force  Armament  Directorate 
Eglln  AFB,  FL  32542-5434 


ABSTRACT 

An  instrumented  package  has  been  accelerated  in  an  electromagnetic  launcher 
to  measure  the  In-bore  acceleration  as  a  function  of  time.  Direct,  continuous 
acceleration  profiles  of  the  launch  show  the  gas  Injection  phase,  electromagnetic 
propulsion,  and  downrange  deceleration  in  a  softcatch  recovery  system.  These 
proof -of -principle  experiments  were  conducted  in  a  50  mm  square  bore  rallgun  and 
utilized  off-the-shelf  components  for  the  In  situ  measurement,  digitization,  and 
storage  of  acceleration  data.  The  1.2  kg  launch  package  was  subject  to  peak 
accelerations  of  nearly  30  kllogees  (2.8  X  10^  m/s^)  in  the  electromagnetic 
propulsion  phase  of  the  launch.  Velocity  and  position  data  obtained  through 
Integration  of  this  data  are  correlated  with  velocity  data  derived  from  con¬ 
ventional  static  diagnostics  (e.g.  ,  magnetic  flux  loops)  to  validate  the  technique. 
The  peak  acceleration  was  slightly  more  than  anticipated  from  the  electrical 
current  delivered  to  the  rallgun,  and  this  deviation  is  examined.  The  presentation 
includes  a  description  of  the  experimental  apparatus,  acquired  data,  a  comparison 
of  the  data  with  code  simulations,  and  suggestions  for  further  application  of 
this  diagnostic  for  in-bore  ballistic  measurements. 
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IN-BORE  ACCELERATION  MEASUREMENTS  WITH 
AN  INSTRUMENTED  RA  LGUN  PROJECTILE 


PRESENTED  AT 

JANNAF  Workshop  on  Electrothermai-Chemical  Modeling  and  Diagnostics 

Aberdeen  Proving  Ground,  MD 
July  9-11,  1991 

PRESENTER 

Donald  M.  Littreil,  USAF,  WL/MNSH 


CONTRIBUTORS 

Keith  A.  Jamison,  SAIC 
Glenn  E.  Rolader,  SAIC 
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PRESENTATION  OUTLINE 


•  PROGRAM  OVERVIEW 

•  EQUIPMENT  r. 

•  PRELIMINARY  TESTS 

•  EML  TESTS 

•  DATA  ANALYSIS 

•  SUMMARY 


JANNAF  ETC  Workshop  July  1 991 


D.M.  Littrell 


PROGRAM  DEFECTION 


In-Bore  Instrumentation/ Diagnostics  (IBID) 


s  The  use  of  self-contained,  on-board  instrumentation  to 
measure  and  explore  the  projectile  environment  within  an 
electromagnetic  gun  during  launch. 


JANNAF  ETC  Workshop  July  1 991 


D.M.  Littrell 


PROGRAM  GOALS 


•  CONSTRUCT  A  PROJECTILE  TO  ACQUIRE  AND  STORE  DATA  DURING 
ACCELERATION  BY  AN  ELECTROMAGNETIC  LAUNCHER 

•  EXPAND  THE  KNOWLEDGE  BASE  ON  PROJECTILE  ACCELERATION  VIA 
THE  APPLICATION  OE  MAGNETIC  FORCES 

•  FURTHER  THE  DEVELOPMENT  OF  GUN-LAUNCHED,  SMART  PROJECTILES 
TO  VERY  HIGH  VELOCITIES 


JANNAF  ETC  Workshop  July  1991 


D.M.  Littrell 
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IBID  CRITICAL  ISSUES 


•  SURVIVABILITY  OF  ELECTRICAL  COMPONENTS  AND  SENSORS  IN  THE 
EML'S  TRANSIENT  ELECTROMAGNETIC.  FIELDS 

•  ELECTRONIC  HARDENING  AGAINST  THE  EML'S  HIGH  ACCELERATION  (G) 
AND  TIME  RATE  OF  CHANGE  OF  ACCELERATION  (G-DOT) 

•  DEVELOPMENT  OF  A  SOFT  CATCH  SYSTEM  TO  STOP  A  HIGH  VELOCITY 
PROJECTILE  WITHOUT  DAMAGE 


JANNAF  ETC  Workshop  July  1991 


D.M.  Linreli 
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PROJECTILE/EML  INTERFACE  ISSUES 


CHARACTERIZATION  OF  IN-BORE  GUN  ENVIRONMENT 

•  JERK  (G-DOT)  DURING  CURRENT  RISE 

•  MAGNETIC  FIELDS 

•  BALLOTING,  EXCESSIVE  LATERAL  GEES 

•  NEGATIVE  G-DOT  IN  FORCING  CONE 

•  NEGATIVE  G-DOT  AT  CROWBAR  AND  EXIT 

•  MUZZLE  WHIP 


jANNAF  ETC  Workshop  JuJv  1991 


D.M.  Littreil 
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LAUNCH  PACKAGE  SCHEMATIC 


JANNAF  ETC  Workshop  July  1991 


D.M.  Littrell 
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INSTRUMENTATION 


ACCELEROMETER  ENDEVCO  MODEL  7270A-60K 

ENDEVCO  MODEL  7270Ar2G0K 

RECORDER/  lES  MODEL  31 

COMPUTER  INTERFACE 

GEE  SWITCH  AGCUDYNE  MODEL  100050-350 

POWER  SUPPLY  DURACELL  MODEL  MN1604  (REPACKAGED) 


JANNAF  ETC  Workshop  July  1 991 


D.M.  Uttrell 
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ACCELEROMETER  SPECIFICATIONS 


Model » 

Range  (kgees) 

Overrange  Limits  (kgees) 
Non-Linearity  &  Hysteresis  (%) 
Transverse  Sensitivity  (%) 
Frequency  Response  (kHz) 

Mounted  Resonant  Frequency  (kHz) 
Mass  (grams) 


7270A-60K 

7270A-200K 

60 

200 

180 

200 

±2 

+2 

5 

5 

DC  to  100 

DC  to  150 

700 

1200 

1.5 

1.5 

JANNAF  BTC  Workshop  July  1 991 


D.M.  Littrell 


HOWITZER  TEST  SETUP 


Hardened  1  Foot 
Concrete  Walls 


155mm  Howitzer 


2x4  Wood 
Target 
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RECORDER  SPECIFICATIONS 


Number  of  Channels  .... 

Bandwidth . 

Input . 

Filter  . . . . 

Amplitude  Resolution  ... 

Stmage  Capacity . . 

Sample  Rate 
Pre-trigger  Storage  ....... 

Triggers  . . . 

Supply  Voltage  ............. 

Supply  Current 

Standby . 

Data  Acquisition 
Dau  Retention  .. 
♦Selectoble  at  Factory 


.  1 

. DC  to  300  KHz* 

. Instrumentation  Amplifier 

. 4  Pole  Butterworth* 

. 8  Bits  (2S6  digital  steps) 

. 32K  Bytes  (32K  x  8)* 

. Up  to  2.5  MHz* 

. .  1762  Bytes 

Switch  Gosure  or  5  Volt  Pulse 
. 8  to  15  Volts 


.  40  micioamps  (typical) 
i!..  60  milliamps  (tj^ical) 
.  40  microamps  (typical) 


JANNAF  ETC  Workshop  July  1991 


D.M.  Littrell 


Acceleration,  Gee's 


HOWITZER  ACCELERATION  DATA 


-63.4K 

-40.2K 

-20.1  K 

0 

20.1  K 

37.0K 

0.0  7.9  15.8  23.6  31.5 

Time,  milliseconds 


1st  concrete 
wall 


Projectile  leaving  the 
gun  barrel  sleeve 


Projectile  acceleration 
in  gun  barrel 


2nd  concrete 
wall 


JANNAF  ETC  Workshop  July  1 991 


D.M.  Littroll 


STATIC  EML  SETUP 

Transient  Field  Effects  Analysis 


Armature 


Current 


JANNAF  ETC  Workshop  July  1991 


D.M.  Littrell 
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2.0 

Volts 

1.0 

0 

-1.0 
-2.0 

1,56  1,80  2.04 

Milliseconds 

”  Perpendicular’ 


JANNAF  ETC  Workshop  July  1 991 


ED  NOISE 

it  Board  Orientation 


Parallel” 


D.M.  Littrell 
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SUMMARY  OF  EML  ACCELEROMETER  TESTS 


Test  Conditions 
Maximum  current  (kA) 
Distance  traveled  (m) 
Injection  velocity  (m/s) 

EML  velocity  (m/s) 

Injector  acceleration  (kgee) 
EML  peak  acceleration  (kgee) 
Projectile  mass  (kg) 

Kinetic  energy  (kJ) 

Bore  size  (mm) 

Armature  type 


PKD 

Boeing 

550 

150 

<1.3 

<0.18 

.  0 

<300 

— 

18 

N/A 

12 

15 

50 

25 

IBID197 

IBID235 

821 

941 

>7 

>7 

149 

162 

264 

360 

3.5 

3.5 

20 

28 

1.17 

1.15 

48 

74 

51 

51 

plasma 

plasma 

JANNAF  ETC  Workshop  Julv  1 991 


O.M.  Littrell 


C8C 

Acceleration  (m/s/s) 


IBID197  ACCELERATION  RECORD 
(a)  Raw  Data  and  (b)  Filtered  Data  (12  kHz) 


400000 
300000 
200000 
100000 
0 

-100000 
-200000 
-300000 
-400000 

0  .005  .01  J015  JOZ  .025  .03  .035  .04  .045  .05  .055 
Time  (seconds) 


G 

o 

w 

(6 

U 

% 

g 

< 


400000 

300000 

200000 

100000 

0 

-100000 

-200000 

-300000 

-400000 


: - 

(b) 

filtered 

r— T  . 1 

• 

- 1 

• 

[ 

_ 

I 

1 

n 

— —7 

HP 

iwfpr 

V 

> 

“""1 

' 

I  i 

0  .005  .01  4)15  .02  .025  .03  4)35  .04  .045  4)5  .055 
Time  (seconds) 


JANNAF  ETC  Workshop  July  1 991 
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CONCLUSIONS 


SUMMARY 

•  TWO  GOOD  DATA  SETS 

•  AGREEMENT  WITH  GAS  GUN  SIMULATION 

•  AGREEMENT  WITH  PROBE  VELOCITY  AND  POSITION  MEASUREMENTS 

•  POTENTIAL  FOR  STUDY  OF  ARMATURE  DYNAMICS 

FUTURE  PLANS 

•  NEAR  TERM  MEASURE  ACCELERATION  AT  CONSTANT  CURRENT 

•  MID  TERM  DEVELOP  MORE  ADVANCED  PAYLOADS  WHICH 

MEASURE  OTHER  IN-BORE  PHENOMENA 

•  FAR  TERM  DEMONSTRATE  A  HIGHLY  INSTRUMENTED  PACKAGE  TO 

SIMULTANEOUSLY  MEASURE  AND  RECORD  MULTIPLE 
CRITICAL  PARAMETERS  RELATING  TO  SMART 
PROJECTILES 
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AGENDA 


8:30 

8:45 


9:00 

9:30 
10:15 
10:35 
10 ;  50 
11:15 

11:40 

12:05 


1:15 
1:30 
1:45 
2:15 
2;  30 

3:00 

3:30 

4:00 


Tuesday,  July  9,  1991 


Registration  ''' 

Welcome  I. 

Administrative  Remarks  G. 

SESSION  I:  SERVICE  USES  OF  ETC  GUNS  &  PLASMAS 
Chairman:  V.  Morelli,  EAPO 


"National  Electric  Gun  Reviews  and  Implications 
for  the  Army's  ETC  Gun  Technology  Program" 

"Electrothermal  Gun  Demonstration  Program" 

"Electrothermal  (ET)  Gun  Program" 

Break 

"Plasma  Discharge  in  the  Electrothermal  Gun" 

"Diagnostics  and  Modeling  of  an  Electrothermal 
Plasma  Source  Experiment  (SIRENS)" 

"Finite  Element  Analysis  of  Engineering 
Electromagnetics  of  ETC  -Gun^" 

Lunch 


V  ¥• 

C. 

s. 

J. 

J. 

0. 

R. 


SESSION  II:  PROPELLANTS 
Chairman:  C.  Daapler»  NSSC 

"Array  Atlernate  ETC  Propellant  Program"  D. 

"Overview  of. Solid  Propellant  ETC  Guns"  A. 

"Overview  of  Gell/Slurry  Propellant"’  A. 

Break 

•  '  ■  •  i 

"Electrothermal -Chemical  (ET-C)  Alternate 

Propellant  Systems  Investigation  and  Study 

Effort"  H. 


"What's  Wrong  with  Therraochemical  Codes 


Applied  to  ETC  Systems?"  .  E. 

"Assessing  ETC  Performance  for 

Systems  Integration"  L. 

Adj  ourn 


May  . 
Wren 


Oberle 

Dampler 

Fowler 

Powell 

Gilllgan 

Hankins 

Boggavarapu 


Downs 

Juhasz 

Bracuti 


McElroy 

Freedman 

Harris 
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8:00 

8:05 

8:35 

9:05 

9:40 

9:55 

10:25 

10:45 

11:15 


12:15 

12:45 

1:15 

1:30 

2:00 

2:30 

3:00 

4:00 


Wednesday.  July  10,  1991 


SESSION  III:  MIXING  &  CONTROL 
Chairmah:  D.  Downs,  ARDEC 

Administrative  Remarks 

"Electrothermal-Chemical  Gun  Program" 

"Diagnostics  Development  for  the  ETC 
Program" 

"Development  of  an  Upvind/Impllcit  Computational 
Model  for  the  Advancement  of  Pitmy  ETC  Guns" 

Break 

"Recent  Advances  in  CAPja  Gtjn  Modeling" 

"30-MM  ETC  Ballistic  Diagnostic  Facility" 

"Numerical  Simulation  of  the  Interior  Ballistic 
Processes  in  an  ETC  Gun" 

Lunch 

SESSION  IV:  MIXING  &  CONTROL 
Chairman:  S.  Vosen,  SMLL 

"Finite-Element  Modeling  of  Electrothermal- 
Chemical  Guns" 

"Special  Diagnostics  and  Instrumentation" 

Break 

"First  Principles  Modeling  of  a  DNA  60iiim  ETC 
Gun  Design" 

"Physics  of  ETC  Plasma-Fluid  Interactions" 

"Observations  and  Modeling  of  Fundamental 
Electrothermal  Gun  Phenomena" 

Adj  ourn 

Bus  leaves  from  Sheraton  Inn  to  the  Inner  Harbor 


G.  Wren 

R.  Woodfln 

D.  Sweeney 

S .  Dash 

D.  Cook 

K.  White 

K.  Kuo 
F.  Cheung 


N.  Wlnsor 
R.  Richardson 

CC.  Hsiao 
B .  Kashiva 

H.  Davis 
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Thursday,  July  11,  1991 


8:00 

8:30 

9:00 

9:30 

9:45 

10:15 

12:00 


SESSION  V:  LESSONS  LEARNED  FROM  OTHER  FIELDS 
Chairman:  J.  Gllligan,  NC  State  U. 

"Electrothermal-Chemical  Gun  Modeling" 

"Rallgun  Research  Relevant  to  Electrothermal  Guns" 

"In-Bore  Position  and  Velocity 
Measurement  Techniques" 

Break 

"In-Bore  Acceleration  Measurements  with  an 
Instrumented  Rallgun  Projectile" 

Group  Discussion  and  Wrap-up 

AdJ  ourn 


D.  King 
.  Batteh 

R.  Bartsch 

.  Llttrell 
.  Wren 
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Eglin  AFB,  FL  32542-5000 

Aberdeen  Proving  Grourrd 


<uncim.  onhr)i  Commander 

U.S.  Army  Armament,  Munitions 
and  Chemical  Command 
ATTN:  AMSMC-IMF-L 
Rock  Island.  IL  61299-5000 

1  Director 

U.S.  Army  Aviation  Research 
and  Technology  Activity 
ATTN:  SAVRT-R  (Library) 

M/S  219-3 

Ames  Research  Center 
Moffett  Field,  CA  94035-1000 


2  Dir,  USAMSAA 
ATTN;  AMXSY-D 

AMXSY-MP,  H.  Cohen 

1  Cdr,  USATECOM 
ATTN:  AMSTE-TC 

3  Cdr,  CRDEC,  AMCCOM 
ATTN:  SMCCR-RSP-A 

SMCCR-MU 

SMCCR-MSI 

1  Dir.  VLAMO 

ATTN:  AMSLC-VL-D 


10  Dir.  BRL 

ATTN:  SLCBR-DD-T 
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Defense  Nuclear  Agency 
ATTN:  John  Farber 
Dr.  fl.  Rohr 
6801  Telegraph  Road 
Alexandria,  VA  22310 

ODDR&E/R&AT 
ATTN:  John  Appel 
R,  Mentz 

Pentagon.  Room  3D1089 
Washington,  DC  20301-3080 

U.S.  Army  Armament  Research,  and 
Development  and  Engineering  Center 
ATTN:  SMCAR-AEE,  Arthur  Bracuti 
Picatinny  Arsenal,  NJ  07606-5000 

U.S.  Army  Armament  Research,  and 
Development  and  Engineering  Center 
ATTN:  SMCAR-AEE-B,  Bldg.  382 
Dr.  Anthony  Beardell 
Donald  Chiu 
Dr.  David  Downs 

Picatinny  Arsenal,  NJ  07806-5000 

U.S.  Army  Armament  Research,  and 
Development  and  Engineering  Center 
ATTN;  SMCAR-AEE’BR,  Bldg.  382 
Dr.  Lee  Harris 

Picatinny  Arsenal,  NJ  07806-5000 

U.S.  Army  Armament  Research,  and 
Development  and  Engineering  Center 
Electric  Armaments  Program  Office  (EAPO) 
ATTN:  SMCAR-FSC,  Bldg.  329  Annex 
Gregory  Ferdinand 
MAJ  R.  Lunberg 
William  Morelli 

Fire  Support  Armaments  Center 
Picatinny  Arsenal,  NJ  07806-5000 

U.S.  Army  Armament  Research,  and 
Development  arxl  Engineering  Center 
ATTN;  SMCAR-FSC.  Hildi  Libby 
Picatinny  Arsenal,  NJ  07806-5000 

Army  Research  Office 
ATTN;  Dr.  David  Mann 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709 


1  Department  of  the  Navy 

Naval  Sea  Systems  Command 
ATTN:  CDR  Craig  Dampier 
CSEA  06  KR12 
Washington,  DC  20362-5101 

1  Naval  Surface  Warfare  Center 
ATTN:  Thomas  Doran 
Code  G33 

Dalgren,  VA  22448-5000 
1  WL/MNSH 

Electromagnetic  Launcher  Technology 
Branch 

ATTN;  Donald  M.  Littrell 
Hypervelocity  Research  Complex 
Wright  Laboratory/Armament  Directorate 
Eglin  Air  Force  Base,  FL  32542-5434 

1  WUMNSH 

ATTN:  Ai  Young 

Eglin  Air  Force  Base,  FL  32542-5434 

3  Los  Alamos  National  Laboratory 
ATTN:  Dr.  R.  Richard  Barlsch 
Dr.  Harold  Davis 
Dr.  James  Trainor 
Group  P'1,  MS  E526 
Los  Alamos,  NM  87545 

1  Los  Alamos  National  Laboratory 
ATTN:  Dr.  Bryan  Kashiwa 
Fluid  Dynamics  Group,  MS  B216 
Los  Alamos,  NM  87545 

1  Sandia  National  Laboratories 

Advanced  Projects  Division,  Org.  9123 
ATTN:  Dr.  David  Benson 
Albuquerque,  NM  87185-5800 

1  Sandia  National  Laboratories 
Division  1811 
ATTN:  Dr.  Paul  Cahill 
Albuquerque.  NM  87185-5800 

1  Sandia  National  Laboratories 
Advanced  Projects  V 
ATTN;  Dr.  Steve  Kempka 
Albuquerque,  NM  87185-5800 
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1  Sandia  National  Laboratories 

Advanced  Projects  Division  V,  Org.  9128 
ATTN:  Dr.  Ronald  L.  Woodfin 
Albuquerque,  NM  87185-5800 

1  Sandia  National  Laboratories 

ATTN:  Dr.  Thomas  Hitchcock 
Albuquerque,  NM  87185-5800 

1  Sandia  National  Laboratories 

Energetic  Materials  Division,  Org.  8357 
Combustion  Research  Facility 
ATTN:  Dr.  Robert  Armstrong 
Livermore,  CA  94551-0969 

1  Sandia  National  Laboratories 

Combustion  Research  Facility,  Org.  8351 
ATTN:  Dr.  Donald  W.  Sweeney 
P.O.  Box  969 
Livermore,  CA  94550 

1  Sandia  National  Laboratories 

Division  8357 
ATTN:  Dr.  Steven  Vosen 
Livermore,  CA  94551-0969 

1  Maxwell  Laboratories,  Inc. 

ATTN:  Ian  R.  McNab 

VP,  Kinetic-Energy  Technologies 
8888  Balboa  Avenue 
San  Diego,  CA  92123 

1  Maxwell  Laboratories,  Inc. 

ATTN:  G.  Mark  Wilkinson 
8888  Balboa  Avenue 
San  Diego,  CA  92123 

1  S-Cubed  Division  of  Maxwell  Labs 

ATTN:  Eduardo  Waisman 
3398  Carmel  Mountain  Road 
San  Diego,  CA  92121 

1  Princeton  Combustion  Research  Laboratory 

ATTN:  Dr.  Neale  Messina 
4275  U.S.  Highway  One  North 
Monmouth  Junction,  NJ  08852 


1  Olin  Rocket  Research 
ATTN:  Dr.  David  Q.  King 
11441  Willows  Road,  N.E. 

Redmond,  WA  98073 

2  The  Pennsylvania  Slate  University 
Propulsion  Engineering  Research  Center 
Mechanical  Engineering  Department 
ATTN;  Professor  F.  B.  Cheung 

Professor  Kenneth  K.  Kuo 
140  Research  Building  "E" 

University  Park,  PA  16802 

3  North  Carolina  State  University 
Department  of  Nuclear  Engineering 
ATTN:  Dr.  Mohomed  Bourham 

Dr.  John  Gilligan 
Dr.  Orlando  Hankins 
Box  7909 

Raleigh,  NC  27695-7909 

1  Stale  University  of  New  York  at  Buffalo 
Department  of  Electrical  Engineering 
ATTN:  Professor  Dr.  W.  J.  Sarjeant 
312  Bonner  ECE-SUNY/AB 
Buffalo.  NY  14260 

1  University  of  Tennessee  Space  Institute 

ATTN:  Dr.  Dennis  R.  Keefer 
Center  for  Laser  Applications,  MS-14 
Tullahoma,  TN  37388-8897 

1  FMC  Corporation 

Naval  Systems  Division  -  M170 
ATTN:  Dr.  David  Cook 
4800  East  River  Road 
Minneapolis,  MN  55421-1498 

3  FMC  Corporation 

Northern  Ordnance  Division  -  Ml  70 
ATTN:  Dr.  Anthony  GiovanettI 
Dr.  Gene  Johnson 
Jahn  Dyvik 

4800  East  River  Road 
Minneapolis,  MN  55421-1498 
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2  General  Dynamics  Land  Systems  Division 
ATTN:  D.  Rao  L.  Boggavarapu 

P.O.  Box  2074 
Warren,  Ml  48090-2074 

1  General  Dynamics  Land  Systems  Division 

ATTN:  Dr.  Melvin  Widner 
MZ  436-21-14 
P.O.  Box  2074 
Warren,  MI  48090-2074 

1  General  Dynamics  Land  Systems  Division 
ATTN:  Dave  Toepler 
MZ  436-21-19 
P.O.  Box  2074 
Warren,  Ml  48090-2074 

1  General  Electric  Defense  Systems 

Department 

ATTN:  Dr.  John  Mandzy 
Bidg.  43-220 
100  Plastics  Avenue 
Pittsfield.  MA  01201 

3  GT-Devices,  Inc. 

ATTN:  Dr.  J.  Robert  Greig 
Dr.  Derek  Tidman 
Dr.  Neils  Winsor 

5705A  General  Washington  Drive 
Alexandria.  VA  22312 

1  Olin  Ordnance 

ATTN:  Hugh  A.  McElroy 
10101  9th  Street  North 
St.  Petersburg,  FL  33716 

2  Science  Applications  international  Corp. 
ATTN:  Dr.  Jad  Batteh 

Lindsey  Thornhill 

1519  Johnson  Ferry  Road,  Suite  300 
Marietta,  GA  30068 

3  Science  Applications  International  Corp. 
ATTN:  Dr.  Sanford  Dash 

Dr.  Ashwin  Hosangadi 
Neeraj  Sinha 

501  Office  Center  Drive,  Suite  420 
Fort  Washington,  PA  19034-3211 


2  Science  Applications  International  Corp. 
ATTN:  C.C.  Hsiao 

Dr.  Gary  Phillips 
10210  Campus  Point  Drive 
San  Diego,  CA  92121 

1  Science  Applicatbns  International  Corp. 
ATTN:  Dr.  Rex  D.  Richardson 
2109  Air  Park  Road,  SE 
Albuquerque,  NM  87106 

1  Teledyne  Brown  Engineering 
ATTN:  Stuart  Fowler 
P.O.  Box  070007,  MS  19 
Huntsville.  AL  35807-7007 

1  Eli  Freedman  and  Associates 
ATTN;  Dr.  Eli  Freedman 
2411  Diana  Road 
Baltimore,  MD  21209 

1  Paul  Gough  Associates,  Inc. 

ATTN;  Dr.  Paul  S.  Gough 
1048  South  Street 
Portsmouth,  NH  03801 

1  Dr.  Keith  Jamison 

P.O.  Box  1261 
Shalimar,  FL  32579 

1  Dr.  Bruce  van  Duesen 

General  Dynamics  Land  Systems  Division 
P.O.  Box  2074 
Warren,  Ml  48090-2074 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it 
publishes.  Your  comments/answers  below  will  aid  us  in  our  efforts. 

1.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of 
interest  for  which  the  report  will  be  used.) _ _ _ _ _ _ 


2.  How,  specifically,  is  the  report  being  used?  (Information  source,  design  data,  procedure, 
source  of  ideas,  etc.)  _ 


3.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or 
dollars  saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please 
elaborate. 


4.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports? 
(Indicate  changes  to  organization,  technical  content,  format,  etc.)  _ 


BRL  Report  Number  BRL-SP-92 _  Division  Symbol 

Check  here  if  desire  to  be  removed  from  distribution  list.  _ 

Check  here  for  address  change.  _ 

Current  address:  Organization  _ _ 

Address  _ 


Department  of  the  Army 

Director 

U.S.  Army  Bailistic  Research  Laboratory 
ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21005-5066 
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BUSINESS  REPLY  MAIL 
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